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Abstract 


The  feasibility  of  using  magnetic  particulate  susceptor  materials  for  induction  heating 
during  bonding  of  polymer  matrix  composite  materials  is  investigated.  If  properly 
designed,  these  systems  should  rapidly  heat  to  the  particulate  material  Curie  temperature 
and  dwell  at  that  temperature  without  feedback  control.  This  performance  is  only 
possible  by  maximizing  hysteresis  heating  while  eliminating  eddy  current  heating  and 
shielding  effects.  Models  of  eddy  current  heating,  hysteresis  heating,  field  penetration 
thickness,  and  heat  transfer  in  a  particulate  film  are  presented.  These  models  are  then 
used  to  predict  heating  behavior  of  particulate  films  containing  metallic  ferromagnetic 
(nickel),  hard  magnetic  (strontium  ferrite),  and  soft  magnetic  (nickel  zinc  ferrite) 
materials.  The  results  show  that  soft  or  semi-hard  insulating  magnetic  materials  are  best 
suited  for  susceptor  particles  for  a  broad  range  of  potential  military  applications.  Polymer 
films  with  around  10%  volume  fraction  of  these  materials  should  be  capable  of  heating 
rates  up  to  100  °C/s  in  moderate  strength  1-10  MHz  induction  fields. 
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1.  Introduction 


1.1  Joining  and  Repair  of  Composite  Materials.  The  goal  of  creating  a  light¬ 
weight,  rapid-insertion,  highly  lethal,  yet  survivable  ground  force  for  the  Army  of  2010  and 
beyond  has  led  to  an  increased  emphasis  on  using  fiber-reinforced  polymer  matrix  composite 
materials  for  both  structural  and  armor  components.  A  major  hurdle  towards  implementation 
of  these  materials,  however,  is  developing  mature  technologies  for  joining  and  repair  of  com¬ 
posite  materials. 

In  order  to  incorporate  composite  components  into  large  assemblies,  they  must  be  joined 
to  many  different  materials,  including  metals,  ceramics,  and  other  polymeric  composites.  The 
simplest  joining  method  is  mechanical  fasteners  [1].  However  this  approach  results  in  poor 
load  transfer,  stress  concentrations,  high  parts  count,  and  increased  weight.  For  these  reasons, 
most  composites  are  joined  using  surface  bonding  techniques  such  as  adhesives  or  thermo¬ 
plastic  fusion  bonding.  Adhesives  often  provide  adequate  bond  strength  and  allow  bonding 
between  dissimilar  materials,  but  typically  require  long  process  times  and  extensive  surface 
preparations  [2].  Recent  work  has  shown  that  adhesive  cure  times  can  be  shortened  without 
bond  strength  reduction  through  the  use  of  elevated  temperature  cure  cycles  [3].  Thermoplastic 
fusion  bonding  joins  materials  through  the  heating,  melting,  and  consolidation  of  a  thermoplas¬ 
tic  interface  [4,  5].  Thermoplastic  composites  can  be  directly  joined  to  other  thermoplastics, 
metals,  or  ceramics  using  this  technique.  Incorporating  a  fusible  thermoplastic  surface  layer 
onto  thermoset  composites,  such  as  through  the  diffusion  enhanced  adhesive  (DEA)  technique 
[6],  also  allows  application  of  thermoplastic  fusion  bonding  to  thermoset  composites.  As  com¬ 
pared  with  conventional  adhesives,  fusion  bonding  typically  offers  shorter  process  times  and 
little  or  no  surface  preparation. 

For  both  accelerated-cure  adhesives  and  fusion  bonding,  localized  heating  is  required  at 
the  bondline.  The  amount  of  heat  generated  must  be  sufficient  to  raise  the  interface  to  100°  - 
200°  C  for  adhesives,  and  150° -400°  C  for  fusion  bonding  of  engineering  thermoplastics.  Higher 
heating  power  allows  for  faster  heating  rates,  reducing  process  times  and  localizing  thermal 
effects  to  the  interface.  Localized  heating  at  the  interface  is  critical  to  successful  bonding  in 
order  to  minimize  energy  requirements  and  prevent  thermal  degradation  of  the  surrounding 
material  [7].  Direct  interface  heating  also  allows  the  joining  of  thick  parts,  which  are  difficult 
to  heat  using  through-thickness  heating  methods. 

A  number  of  interface  heating  techniques  have  been  developed,  including  ultrasonic  weld¬ 
ing  [8],  spin  welding  [9],  and  resistive  heating  [4,  5].  In  ultrasonic  welding,  the  two  materials 
to  be  joined  are  held  in  contact  and  vibrated  at  high  frequencies,  generating  frictional  heat 
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at  their  interface.  This  technique  is  best  suited  to  small  parts  and  simple  geometries.  Spin 
welding  works  on  similar  principles,  except  that  frictional  heating  is  generated  by  relative  ro¬ 
tational  motion  between  the  two  components.  Both  spin  and  ultrasonic  welding  are  forms  of 
thermoplastic  fusion  bonding,  and  are  not  applicable  to  accelerated-cure  adhesives.  In  resis¬ 
tive  heating,  electrical  current  is  passed  through  an  interface  implant  made  of  metal  screen  or 
continuous  carbon  fibers,  which  rapidly  heats  due  to  Joule  heating.  This  approach  has  shown 
great  promise  for  joining  composites  through  both  fusion  bonding  and  elevated-temperature 
adhesives,  with  rapid  heating  rates,  and  high  control  and  repeatability.  However,  there  are  a 
number  of  limitations  to  this  method.  The  need  to  establish  direct  electrical  contact  with  the 
implant  places  geometric  limitations  on  its  applicability.  Power  and  implant  resistance  limita¬ 
tions  require  large  joints  to  be  made  in  many  small  steps.  The  need  for  direct  application  of 
high  electrical  currents  raises  safety  concerns  for  field  applications.  And  the  presence  of  the 
resistive  implant  at  the  interface  often  weakens  the  final  bond  strength. 

1.2  Induction  Heating.  Induction  heating  has  been  proposed  as  a  means  of  interface 
heating  polymer  matrix  composite  materials  [10-13].  In  this  approach,  a  susceptor  material  is 
placed  at  the  bondline  and  exposed  to  an  alternating  magnetic  field.  The  magnetic  field  can 
induce  heating  in  the  susceptor  material  through  two  mechanisms.  If  the  susceptor  material  is 
conductive,  currents  are  induced  in  the  susceptor,  which  then  heats  due  to  resistive  effects.  This 
mode  of  heating  is  called  eddy  current  heating,  and  is  the  only  source  of  heating  in  nonmagnetic 
susceptor  materials  such  as  stainless  steel  and  carbon  fibers.  If  the  susceptor  is  magnetic, 
hysteresis  losses  associated  with  the  magnetization-demagnetization  cycle  causes  additional 
heating.  This  mode  of  heating  is  called  hysteresis  heating. 

A  major  advantage  of  the  induction  method  is  that  it  is  a  noncontact  technique.  Typically 
the  induction  head,  which  generates  the  magnetic  field,  is  placed  on  the  surface  of  the  parts  to  be 
joined.  The  magnetic  field  penetrates  to  the  interface  where  it  causes  the  susceptor  material  to 
heat.  Since  direct  access  to  the  interface  is  not  needed,  induction  heating  offers  great  geometric 
flexibility,  applicability  to  large  bond  areas,  and  improved  safety  as  compared  with  resistive 
heating  methods. 

A  further  benefit  of  the  induction  joining  technique,  when  implemented  with  thermoplastic 
fusion  bonding,  is  that  it  could  improve  the  repairability  of  composite  structures  in  the  field. 
Because  the  susceptor  material  remains  embedded  at  the  interface,  at  any  time  the  interface 
can  be  reheated  through  application  of  a  magnetic  field.  A  joint  which  becomes  damaged  at 
the  bondline  can  be  rejoined  by  remelting  and  reconsolidating  the  interface,  potentially  recov¬ 
ering  the  original  bond  strength.  Components  which  are  damaged  away  from  the  bondline  can 
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be  repaired  by  remelting  the  interface,  removing  the  damaged  component,  and  rebonding  a 
replacement  component. 

1.3  Particulate  Susceptor  Materials.  Traditional  induction  susceptor  materials 
include  metal  screens  and  carbon  fiber  implants.  However,  a  novel  approach  which  offers 
many  advantages  over  traditional  susceptor  materials  is  to  disperse  magnetic  microparticles  at 
the  interface  [14].  The  small  size  of  the  particles  reduces  the  mechanical  degradation  of  bond 
strength  associated  with  traditional  macroscopic  susceptor  materials.  Particulate  susceptors 
also  offer  unique  geometric  possibilities,  such  as  flowable  and  moldable  susceptor  materials. 

The  principal  advantage  of  particulate  susceptor  materials,  however,  is  their  unique  heating 
behavior.  If  the  particle  is  small  enough  or  resistive  enough,  there  is  insufficient  current  path 
for  appreciable  resistive  heating.  In  this  situation  all  of  the  heating  is  due  to  hysteresis  effects. 
A  special  property  of  magnetic  materials  is  the  presence  of  a  Curie  temperature,  which  marks 
the  transition  of  the  material  from  a  ferromagnetic  or  ferrimagnetic  state  to  a  paramagnetic 
state.  A  material  which  is  above  its  Curie  temperature  loses  its  magnetic  properties,  and  can 
no  longer  generate  hysteresis  heating.  Therefore,  if  magnetic  particulate  susceptors  can  be 
designed  which  possess  a  Curie  temperature  and  do  not  heat  due  to  eddy  current  effects,  they 
cannot  be  heated  inductively  above  their  Curie  temperature. 

This  property  allows  the  creation  of  specialized  susceptors  which  “automatically”  maintain 
a  desired  processing  temperature  without  the  need  for  external  control.  As  compared  with 
traditional  induction  bonding  techniques,  this  feature  greatly  improves  process  control  while 
reducing  the  necessary  equipment  sophistication  and  operator  expertise.  This  method  also 
provides  an  opportunity  for  multi-material  cures.  By  layering  different  materials  with  specific 
microparticle  compositions  at  each  interface,  a  single  imposed  magnetic  field  can  be  used  to 
simultaneously  heat  two  interfaces  to  different  dwell  temperatures.  A  thermoplastic  fusion 
bonding  technique  utilizing  induction  heating  of  magnetic  microparticles  therefore  possesses 
tremendous  potential  as  a  simple,  flexible,  noncontact,  highly  controllable  joining  process  for 
composite  materials. 

1.4  Objectives  and  Approach.  In  this  report  we  evaluate  the  feasibility  of  using 
particulate  susceptor  materials  for  induction  bonding  and  repair  of  composite  materials.  In 
addition,  we  will  provide  guidance  as  to  what  materials  and  equipment  are  best  suited  for  this 
application,  and  which  research  thrusts  are  most  crucial  towards  implementation  and  sophisti¬ 
cation  of  this  technique.  The  performance  and  design  goals  for  this  application  are: 
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•  Ensure  automatic  temperature  control  by  choosing  a  system  with  high  hysteresis  heating 
and  low  eddy  current  heating. 

•  Enable  application  to  elevated  temperature  adhesives,  which  require  processing  temper¬ 
atures  from  ~  100°fo200°C,  and  all  engineering  thermoplastics,  which  require  process¬ 
ing  temperatures  from  ~  150°fo400°C. 

•  Keep  equipment  power,  size,  and  cost  within  reasonable  bounds  for  widespread  field 
implementation. 

Because  our  objective  is  only  to  provide  general  recommendations  on  particulate  induction 
heating,  all  calculations  are  performed  to  order-of-magnitude  accuracy.  This  level  of  accuracy 
is  sufficient  to  eliminate  unreasonable  approaches,  compare  vastly  different  approaches,  and 
estimate  feasibility.  The  creation  of  quantitatively  accurate  models  for  design  and  control  of 
specific  induction  processes  will  be  presented  in  future  reports. 

There  has  been  considerable  research  on  the  microscopic  mechanisms  which  cause  hys¬ 
teresis  in  magnetic  materials,  most  notably  domain  wall  motion  [15-19]  and  magnetization 
rotation  [20-22].  However,  these  studies  have  yet  to  provide  a  reliable  means  of  quantitatively 
predicting  macroscopic  magnetic  properties  and  behaviors  from  first  principles.  Therefore,  we 
will  only  address  microscopic  processes  in  order  to  get  a  qualitative  understanding  of  some 
macroscopic  behavior.  The  analytical  modeling  of  heating  processes  will  instead  be  based  on 
reported  macroscopic  magnetic  properties. 

This  report  is  organized  as  follows.  In  section  2,  mathematical  models  of  the  relevant 
process  physics  are  derived.  In  section  3,  these  process  models  are  used  to  determine  the  design 
parameters  for  this  application.  In  section  4,  three  specific  material  systems,  representing  three 
distinct  classes  of  magnetic  materials,  are  evaluated  based  on  the  design  parameters.  Sections  5 
and  6  discuss  and  summarize  the  results  of  this  report,  and  section  7  proposes  future  directions 
for  development  of  this  technology. 

2.  Process  Physics  Models 

In  this  section  we  construct  mathematical  models  of  the  relevant  electromagnetic  and  heat 
transfer  phenomenon  which  are  involved  in  the  proposed  induction  heating  approach. 

2.1  Magnetic  Theory  Preliminaries  and  Definitions.  When  a  magnetic  ma¬ 
terial  is  placed  in  a  magnetic  field,  the  magnetic  field  is  locally  augmented  through  atomistic 
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mechanisms  in  the  magnetic  material.  This  behavior  can  be  shown  graphically  through  an 
H  —  M  magnetization  curve,  figure  1 ,  where  H  is  the  applied  magnetic  field  and  M  is  the 
additional  magnetic  field  created  by  the  material.  M  is  conventionally  referred  to  as  the  mag¬ 
netization  of  the  material.  In  a  vacuum,  the  perfect  nonmagnetic  material,  the  magnetization 
M  has  a  value  of  zero  for  any  applied  field  H.  Once  the  material  reaches  its  saturation  magne¬ 
tization  Ms,  further  increases  in  applied  field  H  do  not  increase  the  magnetization  value. 

The  magnetization  M  is  written  in  terms  of  the  applied  field  H  as 


M  =  (fir  —  l)H,  (1) 

where  p,r  is  a  dimensionless  material  parameter  called  the  relative  magnetic  permeability,  or 
simply  “permeability.”  The  relative  permeability  of  a  vacuum  is  unity.  Equation  1  is  actually  a 
simplification  of  the  general  magnetization  curve,  since  clearly  from  figure  1  the  magnetization 
M  is  not  linearly  related  to  the  applied  field  H.  The  relationship  can  be  written  more  generally 
as 


M(H)  =  (fir(H)  -  1  )H.  (2) 

In  magnetic  material  property  tables,  permeability  is  rarely  tabulated  as  a  function  of  the  ap¬ 
plied  field  H.  Instead,  the  initial  permeability  (at  low  H)  or  some  characteristic  permeability 
averaged  over  a  range  of  H  is  given. 
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Figure  2.  Typical  DC  Hysteresis  Curve. 

In  an  alternating  magnetic  field,  the  magnetization  of  the  material  typically  demonstrates 
hysteretic  behavior  as  shown  in  figure  2.  This  figure  represents  the  saturation  hysteresis  loop, 
where  the  saturation  magnetization  of  the  material  is  reached.  The  magnetization  at  zero  ap¬ 
plied  field  is  called  the  remanance  Mr  and  the  applied  field  necessary  to  reduce  the  magnetiza¬ 
tion  to  zero  is  called  the  coercivity  Hc.  These  values  are  also  frequently  tabulated  in  magnetic 
property  data  tables.  Unless  otherwise  noted,  these  values  are  measured  in  a  steady  (zero  fre¬ 
quency)  magnetic  field. 

We  will  now  formalize  our  language  for  discussing  magnetic  materials  and  fields.  In  this 
report  we  will  use  the  terminology  “magnetic  field”  to  refer  to  B  (in  Teslas,  where  1  T  = 
1  We/m2  =  1  J/A  m2),  which  is  proportional  to  the  sum  of  the  applied  “H”  field  (we  choose  not 
to  name  this  quantity)  and  the  material  magnetization  M,  or 

B  =  fx0(H  +  M) ,  (3) 

where  fi0  is  the  magnetic  permeability  of  free  space,  whose  exact  value  is  defined  as 

H0  =  4tt  x  10-7  H/m.  (4) 

The  units  of  the  magnetic  permeability  are  Henries  per  meter,  which  can  be  rewritten  as  1  H/m 
=  1  N/A2.  The  magnetic  field  B  is  essentially  the  total  magnetic  field,  created  by  the  sum  of 
the  applied  field  H  and  the  material  magnetization  M.  Under  our  convention  the  units  of  B 
differ  from  the  units  of  H  and  M.  (This  difference  in  units  is  not  meant  to  imply  any  special 
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physical  difference  between  the  character  of  B  relative  to  H  and  M.  Rather,  it  exists  mainly 
as  an  artifact  of  historical  convention.)  Equation  3  can  be  rewritten  in  terms  of  the  relative 
magnetic  permeability  as 


B  =  fiofirH.  (5) 

Since  the  magnetization  M,  and  therefore  the  magnetic  field  B,  are  not  linearly  related  to  the 
applied  field  H ,  this  equation  can  be  written  more  generally  as 

B(H)  =  /j,0p,r(H)H.  (6) 

In  all  of  our  equations  we  have  written  H,  M,  and  B  as  scalar  quantities.  They  are  more 
correctly  written  in  vector  form,  where  equation  5  becomes 

Bi  —  Hj .  (7) 

Indicial  notation  is  used  here  with  the  Einstein  summation  convention.  Note  that,  in  general, 
the  relative  permeability  of  a  material  is  a  tensor.  For  the  remainder  of  this  paper  we  will 
assume  isotropic  materials  and  uniform  unidirectional  magnetic  fields,  so  that  the  scalar  forms 
of  the  magnetic  equations  can  be  used. 

Section  2.2  will  require  use  of  the  magnetic  flux  $,  defined  for  some  conductive  loop 
perpendicular  to  a  magnetic  field  as  the  product  of  the  magnetic  field  strength  and  the  loop 
area,  or 


$  =  BAh  (8) 

where  Ai  is  the  area  enclosed  by  the  current  loop.  The  magnetic  flux  is  in  many  ways  a  con¬ 
trived  quantity  that  is  only  meaningful  when  used  within  the  context  of  Faraday’s  law,  which 
will  be  presented  in  equation  1 3.  (A  word  of  caution  to  readers  who  are  new  to  electromagnetic 
field  [EMF]  theory  and  terminology.  Note  that  the  “magnetic  field”  B  is  equal  to  the  “magnetic 
flux”  $  normalized  by  some  area  through  which  it  acts.  This  convention  is  exactly  opposite  to 
the  normal  definition  of  “flux”  in  physics  and  engineering.  If  instead  the  following  terminology 
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was  used 


B  — >  magnetic  flux  [We/m2], 

$  — >  magnetic  flow  [We], 

H  ->  magnetic  field  strength  [A/m], 

then  the  vocabulary  would  be  more  consistent  with  other  branches  of  physics.  Some  authors 
have  endeavored  towards  this  end  by  calling  H  the  magnetic  field  strength,  $  the  magnetic 
flux,  and  B  the  magnetic  flux  density.  However,  the  resulting  language  is  perhaps  even  more 
confusing  than  the  original.) 

2.2  Eddy  Current  Heating.  A  conductive  material  placed  in  an  alternating  magnetic 
field  has  alternating  currents  induced  within  it.  These  currents  generate  heat  in  the  material  due 
to  Joule  (resistive)  effects,  a  phenomenon  known  as  eddy  current  heating. 

Consider  a  spherical  particle  of  radius  R0  and  resistivity  pa  in  the  presence  of  a  uniform 
periodic  magnetic  field  of  strength 


B  =  Bp  sino>/, 


(9) 


where  Bp  is  the  field  amplitude  and  u  is  the  frequency  of  the  field  (in  radians  per  second). 
By  assuming  a  uniform  field  and  a  spherical  particle,  the  exact  direction  of  the  field  is  not 
important.  However,  for  calculation  purposes  we  will  define  a  cylindrical  coordinate  system 
such  that  the  magnetic  field  acts  in  the  —  z  direction. 

Because  induction  heating  is  a  volumetric  phenomenon,  we  will  be  performing  a  volumetric 
integration.  The  discretization  of  the  geometry  is  shown  in  Figure  3,  for  a  slice  of  the  sphere 
in  the  r  —  9  plane.  We  will  use  the  notation  R(z)  to  represent  the  radius  of  the  slice,  where 

R(z)  =  (10) 


The  resistance  Rq  of  the  circumferential  strip  of  width  dr  and  height  dz  is 


Rci  =  Pci 


length 


cross-sectional  area 


=  pci 


2nr 

drdz 


(11) 
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Figure  3.  Discretization  of  Sphere  for  Eddy  Current  Model. 


The  magnetic  flux  $  for  this  strip  is  defined  as  the  product  of  the  magnetic  field  strength  and 
the  area  enclosed  by  the  strip 

$  =  (Bp  sintuf)  7rr2.  (12) 

Faraday’s  Law  states  that  the  EMF  generated  in  a  current  loop  is  proportional  to  the  rate  of 
change  of  the  magnetic  flux,  or 


For  our  periodic  magnetic  flux,  the  induced  EMF  is  also  periodic 

V  =  —  nr2u>  (Bp  cos  u>t) . 


(13) 


(14) 


Note  that  the  induced  EMF  is  exactly  out  of  phase  with  the  magnetic  flux,  and  its  highest  values 
occur  when  the  magnetic  flux  is  zero  but  the  rate  of  change  of  the  magnetic  flux  is  maximized. 


The  power  consumed  due  to  Joule  heating  in  the  circumferential  strip  is 


V2 

dP  =  — , 

n  ? 

(7rr2Lu  [Bp  coswf])2 

=  .  2ttp  ’ 

PQ  drdz 

7Tr3LU2B 2  cos2  Ujt 

= - £ - drdz. 

2  pn 

Taking  the  time  average  of  this  instantaneous  power  function  gives 


(15) 

(16) 
(17) 


irr3u2B2 

dP  — - -drdz. 

4pn 


The  total  power  consumed  over  the  entire  strip  is 


(18) 


P  =  j  dP, 

(19) 

[Ro  irr3u2Bl  cos2  ut 

=  2  p  drdz , 

Jo  Jo 

(20) 

7TU )2Bl  [Ro  [R (z)  , 

=  /  /  r3dreh, 

Jo  Jo 

(21) 

nuPB2  fR°  1 
=  R(zfdz. 

Jo  4 

(22) 

Using  equation  10, 


ITU) 

P~  - *- 


2  B2  'R 


f  °(Rl  -  z2fdz , 
Jo 


8pn 

8  p5\ 

ITU}2 


'2B2 


L  t?5 

15  pa  °' 

Normalizing  this  quantity  by  the  volume  of  the  sphere  gives 

p  *<sb;ri 

20  pn 


(23) 

(24) 

(25) 


(26) 
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This  formula  can  be  re-expressed  in  terms  of  the  frequency  /  =  u;/27r  (in  cycles  per  second) 
and  sphere  diameter  D0  as 


Pe  = 


(7 r/ffp£0)2 

20pn  ’ 


(27) 


where  we  have  added  the  subscript  e  to  distinguish  this  quantity  as  the  volume-averaged  eddy 
current  power. 

A  number  of  important  trends  can  be  observed  from  equation  27.  The  eddy  current  power 
increases  with  frequency  squared  and  sphere  diameter  squared,  while  it  decreases  with  increas¬ 
ing  resistivity.  Therefore  we  know  that  the  eddy  current  power  can  be  minimized  by  using 
lower  frequencies,  smaller  particles,  and  higher  resistivity  materials. 


2.3  Hysteresis  Heating.  In  the  presence  of  a  magnetic  field,  a  magnetic  material  un¬ 
dergoes  a  process  of  microstructural  alignment.  The  mechanisms  involved  in  this  realignment 
include  magnetization  rotation  and  domain  wall  motion.  In  most  polycrystalline  materials, 
domain  wall  motion  is  the  “easier”  process,  and  is  responsible  for  most  of  the  magnetization 
changes  within  the  material.  Domain  wall  motion  is  a  highly  irreversible  process,  so  that  sig¬ 
nificant  energy  losses  are  associated  with  material  magnetization.  In  an  alternating  magnetic 
field,  these  losses  occur  continuously  over  time,  resulting  in  power  consumption  which  takes 
the  form  of  heating  in  the  material. 

Figure  4  shows  a  typical  H  —  M  hysteresis  curve  for  a  magnetic  material.  By  gradually 
imposing  a  field  Hp,  reversing  the  field  to  —Hp,  and  then  reimposing  the  field  Hp,  the  mag¬ 
netization  M  follows  this  hysteresis  curve.  The  exact  shape  of  the  curve  is  a  function  of  the 
susceptor  material,  the  susceptor  material  temperature,  and  the  frequency  of  the  applied  field. 
If  the  potential  energy  of  the  system  is  defined  as  fi0MH,  then  the  total  energy  lost  during  one 
complete  hysteresis  cycle,  per  unit  volume  of  material,  is 

Eh  =  no  <(  HdM.  (28) 


If  an  alternating  current  is  applied  to  the  material,  then  the  hysteresis  loop  is  traversed  once  per 
cycle,  so  that  the  power  consumed  is 


Ph  =  fEh. 


(29) 
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Figure  4.  General  Hysteresis  Curve. 

This  equation  is  exact,  provided  the  hysteresis  energy  density  is  known  exactly  for  the  process¬ 
ing  conditions  (temperature,  imposed  field  amplitude,  frequency). 

For  a  specific  material  system,  the  best  way  to  implement  equation  29  is  to  measure  the 
hysteresis  curve  for  the  material  at  the  desired  frequency,  amplitude,  and  temperature.  How¬ 
ever,  it  is  rare  to  find  the  exact  desired  hysteresis  curve  for  a  given  material  in  the  literature. 
Therefore  methods  are  needed  for  estimating  a  specific  hysteresis  curve  area  based  on  available 
information. 

DC  saturation  curves  for  the  three  specific  materials  under  consideration  for  susceptor  par¬ 
ticles  are  given  in  section  4.  No  AC  hysteresis  curves  for  these  materials  could  be  found 
in  the  literature.  Although  in  many  materials  the  hysteresis  energy  density  is  frequency- 
dependent  [23, 24],  the  best  estimate  we  can  make  with  the  available  information  is  frequency- 
independent  hysteresis  energy  density  equal  to  the  DC  hysteresis  energy  density. 

Saturation  hysteresis  curves  are  directly  useful  if  the  applied  field  necessary  to  achieve 
saturation  is  less  than  or  equal  to  the  actual  field  strength  which  can  be  applied  with  feasible 
induction  manufacturing  and  repair  equipment.  This  situation  is  common  with  soft  magnetic 
materials.  If  the  applied  field  necessary  to  achieve  saturation  is  much  greater  than  the  feasible 
equipment  field  strength,  as  is  the  case  with  hard  magnetic  materials,  then  magnetic  saturation 
in  the  material  is  never  reached.  In  this  case  it  is  necessary  to  approximate  a  minor  hysteresis 
loop  for  the  material  based  on  its  major  hysteresis  loop. 
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For  a  saturation  hysteresis  loop,  we  can  write  the  saturation  hysteresis  energy  density  as 

Ehs  =  {3»0MSHS,  (30) 

where  Hs  is  the  field  strength  required  to  mostly  saturate  the  material  and  /?  is  a  constant 
bounded  by  [0,1].  It  is  reasonable  to  assume  that  minor  hysteresis  loops  are  shaped  similarly  to 
the  major  loops,  only  scaled  according  to  the  applied  field  strength.  This  assumption  has  two 
consequences.  First,  j3  is  constant,  so  any  minor  hysteresis  loop  can  be  written  as 


Eh  =  Pn0MpHp 


(31) 


Secondly,  the  peak  magnetization  decreases  linearly  with  the  applied  field  strength,  or 


Mp  Hv 

ms  ~  h: 


(32) 


Combining  equations  30  and  31  yields 


tp  —  p  Mz  Mz 

h  hs  Ms  Hs 


(33) 


which,  when  combined  with  equation  32,  gives 


Eh  =  Ehs 


(HA2 
\HS)  ' 


(34) 


For  all  field  strengths  greater  than  Hs,  Eh  will  be  equal  to  Ehs.  Therefore  we  can  write  equa¬ 
tion  34  more  generally  as 


Eh 


[^(f)2  < tHp<a„ 


Ehs 


if Hp  >  Hs. 


(35) 


For  our  derivations  in  section  3.1,  we  will  find  it  useful  to  write  an  alternate  expression  for 
the  hysteresis  energy  density.  We  define  a  hysteresis  permeability  j2h  as 


fJ'h  =  1  + 


Ehs 

PoH* 


(36) 


13 


Figure  5.  Approximation  of  the  Initial  Permeability  from  a  Saturation  Hysteresis  Curve, 
so  that 


Ehs  =  Hoifih  -  1  )H2S.  (37) 

jlh  is  essentially  the  proportionality  constant  between  Ehs  and  H%,  and  can  be  calculated  from 
a  measured  saturation  hysteresis  loop.  We  have  chosen  to  express  this  hysteresis  constant  in  a 
form  analogous  to  a  relative  magnetic  permeability,  because  it  simplifies  later  derivations  and 
has  a  value  of  unity  for  non-lossy  materials.  Equation  34  can  be  rewritten  as 


Eh  -  -  1  )Hp.  (38) 

Therefore,  the  complete  hysteresis  energy  model  of  equation  35  can  be  implemented  for  a 
material  if  Jlh  and  Hs  are  known. 

We  will  also  find  it  necessary  to  estimate  the  true  magnetic  permeability  of  a  material 
from  its  saturation  hysteresis  loop.  This  parameter  will  be  used  in  derivations  which  need  to 
estimate  the  magnetic  field  strength  produced  by  a  material  for  a  given  applied  field  strength. 
Since  permeability  is  actually  a  function  of  applied  field  strength  (equation  6),  we  need  some 
sort  of  average  permeability  value.  Furthermore,  we  will  need  to  estimate  this  value  based  on 
a  given  saturation  hysteresis  curve.  The  approach  we  will  use  is  to  graphically  estimate  the 
average  permeability,  as  shown  in  Figure  5. 
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2.4  Field  Penetration  Depth.  Induction  heating  is  a  volumetric  effect.  The  eddy 
current  analysis  in  section  2.2  assumed  a  uniform  magnetic  field  throughout  the  particle.  How¬ 
ever,  the  induced  currents  in  the  particle  create  a  magnetic  field  which  opposes  the  applied 
magnetic  field,  a  phenomenon  known  as  the  shielding  effect.  Typically,  the  net  magnetic  field 
is  highest  at  the  surface  of  the  susceptor  and  decreases  exponentially  with  depth.  The  pen¬ 
etration  depth  is  usually  expressed  in  terms  of  a  skin  depth,  8,  which  can  be  calculated  by 
[25] 


(39) 


The  skin  depth  is  defined  as  the  distance  at  which  the  magnetic  field  has  been  reduced  to  around 
37%  (1/e)  of  its  surface  value.  If  the  skin  thickness  is  less  than  the  particle  size,  then  shielding 
effects  are  important. 


2.5  Induction  Head.  Induction  heating  equipment  typically  consists  of  a  power  sup¬ 
ply,  an  induction  head,  and  a  cooling  system.  The  power  supply  sends  current  through  the 
induction  head,  which  is  made  of  a  winding  of  low-resistivity  material  such  as  copper.  The 
shape  of  the  head,  its  position  relative  to  the  susceptor,  and  the  current  flowing  through  the 
head  determine  ihe  magnetic  field  applied  to  the  susceptor.  Power  is  consumed  through  heat¬ 
ing  of  the  susceptor,  due  to  eddy  and  hysteresis  losses,  and  through  resistive  heating  in  the 
head.  The  resistance  of  the  head  must  be  minimized  to  reduce  head  losses.  Because  of  the  high 
currents  necessary  to  induce  significant  magnetic  fields,  the  induction  head  is  often  cooled 
through  water  or  air  cooling.  The  most  common  approach  is  to  make  the  induction  head  out  of 
copper  tubing,  and  continually  circulate  cooling  water  through  the  head. 

In  order  to  model  the  power  loss  in  the  head,  the  resistance  of  the  head  must  be  calculated 
through 

Rq  =  Pn~r^-,  (40) 

Aih 

where  Lih  and  A,h  are  the  length  and  cross-sectional  area  of  the  induction  head,  respectively. 
The  power  loss  in  the  head,  P  =  I2 Rq,  is  then 


For  a  given  current  and  head  resistance,  this  formula  allows  calculation  of  the  power  consumed 
by  the  head.  This  value  can  be  used  to  determine  the  necessary  cooling  system  capacity. 

We  must  also  calculate  the  magnetic  field  produced  by  a  typical  induction  head.  The  mag¬ 
netic  field  produced  by  a  single  current  loop  of  radius  R  at  a  distance  2  along  the  centerline  of 
the  loop  is  [26] 


ii0l  R 2 

2  (R2  +  z2) 3/2  ■ 


(42) 


A  typical  “pancake-coil”  induction  head  consists  of  a  set  of  concentric  current  loops.  The 
magnetic  field  produced  by  this  head  can  be  calculated  by  summing  the  contributions  of  each 
loop  as  calculated  by  equation  42. 


2.6  Heat  Transfer  in  Particulate  Film.  We  will  use  the  terminology  “particulate 
film”  to  refer  to  both  thermoplastic  films  and  adhesives  loaded  with  volumetrically-heating 
particles  because,  neglecting  heats  of  fusion  and  reactive  exotherms,  the  heat  transfer  models 
are  identical  for  both  cases.  These  secondary  heat  transfer  effects  are  expected  to  be  insignif¬ 
icant  for  this  order-of-magnitude  feasibility  study,  but  need  to  be  incorporated  into  any  future 
quantitative  modeling  efforts. 

Since  the  particles  and  the  distance  between  particles  are  small,  we  will  use  a  very  simple 
lumped  analysis  for  modeling  the  heat  transfer  in  the  particulate  film.  For  a  film  consisting  of  a 
volume  fraction  vj  of  particles  of  density  p  and  specific  heat  cp  dispersed  in  a  matrix  of  density 
p  and  specific  heat  cp,  the  power  required  to  heat  the  film  over  a  temperature  change  AT  in  a 
time  At  is 


Pf  =  (vJPcp  +  (1  -  vj)pcp) ,  (43) 

where  Pj  is  the  power  consumed  in  the  film,  per  volume  of  film. 

This  heat  transfer  model  assumes  that  the  heat  which  is  generated  within  the  particles  is 
immediately  and  uniformly  distributed  over  the  entire  film.  This  assumption  is  only  valid  if  the 
spacing  between  the  particles  is  sufficiently  small  such  that  heat  can  be  transferred  across  that 
distance  in  a  time  much  shorter  than  the  process  time.  We  can  calculate  the  critical  spacing 
between  the  particles  by  calculating  the  characteristic  conduction  length  scale  for  the  film,  Lc 

Lc  =  \/  <ycAt ,  (44) 
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Figure  6.  Particle  Dispersion  Geometry  for  Characteristic  Conduction  Length  Calculation. 


where  ac  is  the  thermal  conductivity  of  the  film.  For  a  body-centered  cubic  arrangement  of 
spherical  particles,  Figure  6,  the  particle  volume  fraction  can  be  expressed  as 


vf  = 


Lz  ’ 


(45) 


where  L  is  the  particle  spacing  along  the  primitive  cell  edges.  Since  heat  is  only  generated 
within  the  particles,  the  characteristic  conduction  length  is  roughly  the  distance  L  between 
particles.  Equating  L  =  Lc,  combining  equations  44  and  45,  and  rearranging  yields 


(46) 


Taking  3 /n  ~  1  gives 


Dc  =  v)/Zy/^Kt.  (47) 

For  a  given  particle  volume  fraction  vj  and  film  thermal  diffusivity  ac,  the  particle  size  must 
be  less  than  Dc  in  order  to  keep  the  characteristic  conduction  times  less  than  At. 

3.  Design  Parameters 

Now  that  the  relevant  physics  for  this  problem  have  been  modeled,  design  criteria  can  be 
derived  which  place  bounds  on  acceptable  material  properties  and  equipment  capabilities  for 
processing. 


3.1  Eddy  Heating  to  Hysteresis  Heating  Ratio.  For  automatic  temperature 
control  through  Curie  temperature  effects,  the  eddy  heating  in  the  particles  must  be  signifi¬ 
cantly  less  than  the  hysteresis  heating  in  the  particles.  Using  equations  27  and  29,  this  ratio 
is 


,  _  Pe  ...  (~BpD0Yf 
p  Ph  20 paEh  • 


(48) 


This  relation  for  the  heating  ratio  Rp  is  its  most  accurate  form,  and  should  be  used  whenever 
possible. 

Since  Eh  is  dependent  on  Bp,  equation  48  is  not  in  its  most  primitive  form.  Using  equa¬ 
tions  5  and  38  we  can  write 


B[  (wJ,)1  ft 

Eh  1)*,2’ 


Using  this  relationship  to  eliminate  Bp  and  Eh  from  equation  48  yields 


Rp 


7 r 


DU  Pi 


-Po~ 


20  pa  -  1 

Taking  7r2/20  ~  1  the  power  ratio  can  be  written  as 


Bp  —  PqB  { 


2  ,  ft  1  ~ 

hPh  —  1  Pa, 


f 


or 


(49) 


(50) 


(51) 


rp  =  p0D2or1f, 


where 

nr  _  ft  1 

1 1  —  - 7  —  • 

Ph-lpa 


(52) 


(53) 


T!  is  the  relevant  material  parameter  for  determining  the  ratio  of  eddy  heating  to  hysteresis 
heating.  Note  that  the  power  ratio  decreases  (less  eddy  heating)  for  materials  with  higher 
resistivities,  smaller  particles,  and  lower  frequencies.  The  power  ratio  becomes  infinite  for  a 
hysteresis  permeability  of  1.  This  case  corresponds  to  a  material  with  no  hysteresis  losses. 
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Figure  7.  Power  Ratio  Plot. 

We  can  use  equation  52  to  create  a  guideline  for  designing  an  induction  heating  process 
with  negligible  eddy  heating.  We  will  impose  the  (somewhat  arbitrary)  condition  that  eddy 
heating  is  insignificant  if  Rp  <  0.01.  We  can  then  plot  D0  as  a  function  of  Ti  and  /  for  this 
value  of  Rp.  The  resulting  plot.  Figure  7,  shows  the  limiting  particle  sizes  for  insignificant 
eddy  heating  as  a  function  of  the  material  and  imposed  frequency. 

3.2  Dimensionless  Penetration  Depth.  In  section  2.4  we  presented  the  relation¬ 
ship  for  the  magnetic  field  penetration  depth  as  a  function  of  the  material  parameters  and  ap¬ 
plied  field.  If  the  penetration  depth  is  comparable  to  or  less  than  the  particle  dimension,  then 
shielding  effects  are  significant.  In  this  case  both  our  eddy  and  hysteresis  models  will  be  in¬ 
accurate  and  can  no  longer  be  applied,  since  they  both  assumed  a  uniform  applied  field.  The 
actual  field  inside  the  particle  will  be  reduced  due  to  shielding  effects,  and  heating  will  be  di¬ 
minished.  Therefore  we  wish  to  keep  the  particle  radius  less  than  the  penetration  thickness. 
Defining  a  dimensionless  penetration  thickness  ratio  Rs  =  28/ D0,  and  using  equation  39,  we 
can  write 


R .  = 

D0  V 


(54) 
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Figure  8.  Dimensionless  Penetration  Thickness  Plot. 


If  the  dimensionless  penetration  thickness  is  large  then  shielding  effects  are  insignificant.  Since 
both  heating  power  and  shielding  effects  increase  as  frequency  increases,  shielding  mecha¬ 
nisms  will  probably  provide  an  upper  bound  on  the  feasible  process  frequency. 

Taking  2/y/n  ~  1,  equation  54  can  be  rewritten  as 


(55) 


where 

r2  =  x[¥.  (56) 

V  /*r 

T2  is  the  relevant  material  parameter  for  determining  the  dimensionless  penetration  depth.  We 
can  use  this  equation  to  create  a  guideline  for  designing  an  induction  heating  process  with 
negligible  shielding  effects.  We  will  impose  the  (somewhat  arbitrary)  condition  that  shielding 
is  insignificant  if  Rs  >  10.  We  can  then  plot  D0  as  a  function  of  T2,  and  /  for  this  value  of 
Rs.  The  resulting  plot,  Figure  8,  shows  the  limiting  particle  sizes  for  insignificant  shielding  as 
a  function  of  the  material  and  imposed  frequency. 
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Figure  9.  Characteristic  Heating  Rates. 


3.3  Heating  Rates.  The  time  necessary  to  reach  the  processing  temperature  must  be 
small  enough  to  allow  for  rapid  bonding  and  repair.  We  can  rewrite  equation  43  as 

dT 

—  =  VfPh  (vfpcp  +  (1  -  v/)pcp)_1 .  (57) 


The  term  dT/dt  is  the  characteristic  heating  rate  and  vjPh  is  the  particle  heating  power  due 
to  hysteresis  effects  only,  assuming  negligible  eddy  heating.  The  extra  factor  of  vj  is  included 
because  the  hysteresis  heating  is  generated  only  in  the  film  volume  occupied  by  the  particles. 
Combining  this  equation  with  equation  29  and  rearranging  yields 


dT  rEh 

__ _  ~  -f  n 

dt  pCp 


1  + 


(58) 


This  equation  allows  us  to  calculate  the  characteristic  heating  rate  for  a  given  material,  partic¬ 
ulate  volume  fraction,  and  applied  field  frequency.  Figure  9  shows  characteristic  heating  rates 
for  different  heating  times  and  temperature  rises.  The  range  of  100°  —  400°C  encompasses 
the  expected  range  of  processing  temperature  changes  for  elevated-temperature  adhesives  and 
engineering  thermoplastics.  A  heating  time  of  100  s  is  considered  reasonable.  Therefore  we 
need  to  achieve  heating  rates  of  at  least  1°  -  10°C/s. 
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A  number  of  important  trends  are  revealed  by  equation  58.  The  heating  rate  is  proportional 
to  the  field  frequency.  Therefore  increasing  the  frequency  by  an  order  of  magnitude  should 
also  increase  the  heating  rate  by  an  order  of  magnitude.  Increasing  the  hysteresis  loop  energy 
density  also  proportionally  increases  the  heating  rates.  However,  it  is  very  easy  to  increase  the 
field  frequency  by  orders  of  magnitude  through  equipment  selection  and  control,  while  there  is 
no  systemic  means  of  increasing  hysteresis  energy  density. 

The  effect  of  particle  volume  fraction  on  heating  rates  can  be  investigated  by  writing  the 
ratio  of  the  heating  rate  for  general  vj  to  the  particle  heating  rate  (v j  =  1.0)  as 

( dT/dt)V{  _  1 

m  ”  l+il/vj-WCr/pCrY  (59> 

This  equation  gives  the  dimensionless  heating  rate  as  a  function  of  the  particle  volume  fraction 
and  the  specific  heat  ratio  pcpjpcv.  Figure  1 0  shows  the  dimensionless  heating  rate  as  a  function 
of  volume  fraction  for  three  different  specific  heat  ratios.  Note  that  when  the  matrix  specific 
heat  is  less  than  the  particle  specific  heat  (pcp/ pcp  =  0.1),  which  is  the  case  most  common  for 
magnetic  particles  in  polymer  matrices,  significant  changes  in  heating  rate  with  volume  fraction 
are  noted  at  low  particle  volume  fractions.  However,  at  higher  volume  fractions,  little  change 
in  heating  rate  is  noted  as  particle  volume  fraction  is  increased.  This  behavior  is  advantageous, 
allowing  significant  film  heating  with  relatively  low  particle  volume  fractions. 

3.4  Curie  Temperature.  For  systems  with  negligible  eddy  heating,  the  Curie  tem¬ 
perature  represents  the  dwell  temperature  of  the  material  during  induction  heating,  assuming 
the  inductive  heating  power  is  sufficient  to  overcome  thermal  losses  and  heat  the  susceptor  to 
this  temperature.  Since  the  processing  temperature  of  most  elevated-temperature  adhesives  and 
engineering  thermoplastics  is  between  100°  and  400°C,  only  materials  with  Curie  temperatures 
in  this  range  should  be  considered  for  the  particulates. 


4.  Calculations  for  Specific  Material  Systems 

In  this  section  we  perform  quantitative  calculations  of  the  design  variables  for  specific 
equipment  and  material  configurations.  In  section  4.1,  equipment  capabilities  are  quantified 
and  some  preliminary  heat  transfer  calculations  are  performed.  Sections  4.2-4.4  then  pro¬ 
vide  characteristic  heating  behaviors  for  three  different  materials:  nickel,  strontium  ferrite,  and 
nickel  zinc  ferrite,  representing  elemental  ferromagnetic,  hard  ferrite,  and  soft  ferrite  magnetic 
materials,  respectively. 
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Figure  10.  Dimensionless  Heating  Rates  as  a  Function  of  Particle  Volume  Fraction  v j  and 
Matrix/Particle  Specific  Heat  Ratio  pcv/  pcv. 

4.1  General  Equipment  and  Matrix  Properties.  Although  many  induction  coil 
geometries  are  possible,  for  our  benchmarking  purposes  we  will  consider  a  standard  5-tum  pan¬ 
cake  coil  made  out  of  copper  tubing  with  1 0  mm  outside  diameter  and  2  mm  wall  thickness.  We 
will  approximate  the  coil  as  5  concentric  circles  of  radius  2,  4,  6,  8,  and  10  cm.  The  effective 
length  and  cross  sectional  area  of  this  head  are 

Lih  =  27t(20  +  40  +  60  +  80  +  100)  mm  Aih  —  ^-7r(102  —  62)  mm2.  (60) 

A  typical  induction  unit  passes  100  A  of  current.  Using  a  copper  resistivity  value  of  1.7  x 
10~8  Cl  m  [27]  and  the  resistive  heating  power  equation,  equation  41,  the  head  consumes  a 
total  power  of  P^  =  6.4  W.  Since  we  expect  to  use  power  supplies  with  kW  capacity,  the 
power  consumed  by  the  head  is  negligible. 

For  our  pancake  coil  geometry,  and  assuming  the  coil  is  located  2  cm  above  the  susceptor 
material,  equation  42  predicts  a  field  strength  of 

B  =  4.44  x  10~5/,  (61) 
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Table  1.  Typical  Polymer  Matrix  Thermophysical  Properties. 


Property 

Value 

Density  (p) 

1000  kg/m3 

Specific  heat  (cp) 

lOOOJ/kgK 

Thermal  conductivity  ( k ) 

0.2  W/m  K 

where  /  is  in  Amps,  and  B  is  in  Teslas.  For  a  current  of  100  A,  the  magnetic  field  is  B  =  4.44  x 
10~3  T.  Converting  this  value  to  an  H  value  using  H  =  B/p0,  the  field  value  is  3500  A/m. 
This  value  will  be  used  as  our  maximum  field  strength  in  all  subsequent  calculations. 

Commercial  induction  heating  equipment  usually  operates  in  the  range  of  1  kHz  to  1 0  MHz. 
However,  higher  frequency  equipment  can  be  manufactured  if  needed. 

For  our  film  heat  transfer  model,  the  thermal  properties  of  the  suspending  medium  are 
needed.  These  properties,  typical  for  adhesives  and  engineering  thermoplastics,  are  given  in 
Table  1 .  Using  these  matrix  properties  as  the  overall  film  properties  (a  conservative  estimate), 
and  setting  the  characteristic  conduction  time  as  1  s,  equation  47  gives  the  critical  particle 
diameter  for  uniform  heating  as 

Dc  =  v)/3J ^  =  (4.47  x  10->}/3  m.  (62) 

V  PCp 

This  critical  particle  diameter  is  plotted  as  a  function  of  volume  fraction  in  Figure  1 1 .  Note  that 
even  at  particle  volume  fractions  of  only  1%,  the  critical  particle  diameter  is  ~  100  pm.  Ac¬ 
cording  to  this  calculation,  if  the  particles  are  less  than  100  in  diameter,  heating  throughout 
the  film  should  be  uniform.  At  higher  volume  fractions  this  critical  particle  size  increases  to 
~  300  pm.  Since  this  calculation  is  only  guaranteed  to  be  accurate  within  one  or  two  orders  of 
magnitude,  it  is  probably  best  to  use  particles  no  larger  than  1-10  pm. 

4.2  Nickel  (Elemental  Ferromagnetic).  Of  the  five  elemental  ferromagnetic  ma¬ 
terials,  iron,  cobalt,  nickel,  gadolinium,  and  dysprosium,  nickel  is  the  only  one  whose  Curie 
temperature,  358° C,  is  within  the  bounds  of  thermoplastic  processing.  As  a  powder  it  is  also 
readily  available  and  inexpensive.  Table  2  lists  some  of  the  relevant  magnetic  and  thermophys¬ 
ical  properties  for  nickel.  Figure  12  shows  a  DC  saturation  hysteresis  loop  for  3  y,m  nickel 
powder,  measured  by  Suwanwatana  [28].  The  magnetic  properties  in  Table  2  were  calculated 
from  this  DC  hysteresis  curve.  Ehs  is  the  area  of  the  loop,  Hs  is  applied  field  strength  necessary 


24 


Figure  11.  Critical  Particle  Diameter  to  Maintain  Thermal  Conduction  Times  on  Order  of  1  s, 
as  a  Function  of  Volume  Fraction. 

to  mostly  saturate  the  material,  and  fir  is  measured  graphically,  as  suggested  in  section  2.3.  The 
fir  value  of  8.0  suggests  that  at  our  expected  maximum  applied  field  strength,  3500  A/m,  the 
material  magnetization  will  be  Bp  =  0.03  T. 

Using  equation  56  and  the  values  of  pa  and  fir  from  Table  2,  the  penetration  thickness 
material  parameter  T2  has  a  value  for  nickel  of 

T2  =  9.2  x  10-5  01/2m1/2.  (63) 

From  Figure  8,  we  can  see  that  for  this  value  of  X2  skin  effects  may  become  important.  For 
example,  at  1  MHz,  the  particles  must  be  smaller  than  10  pm  for  shielding  to  be  unimportant. 
At  100  MHz,  the  particles  must  be  smaller  than  1  pm. 

Using  equation  36  and  the  values  of  Ehs  and  Hs  from  Table  2,  we  calculate  fih  =  2.0. 
Substituting  this  value  with  the  values  of  pa  and  fir  into  equation  53,  the  heating  ratio  parameter 
T  i  has  a  value  for  nickel  of 


Ti  =  9.4  x  108  H^m"1. 


(64) 
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Table  2.  Nickel  Thermophysical  and  Magnetic  Properties.  All  values  above  horizontal  line 
from  Ozisik  [29]  except  electrical  conductivity,  from  Weast  [27].  All  values  below  line  deter¬ 
mined  graphically  from  Figure  12. 


Property  Value 

Density  ( p )  8900  kg/m3 

Specific  heat  (cp)  430  J/kg  K 

Electrical  resistivity  (pa)  6.8  x  10-8  f)  m 


Curie  temperature  358°C 

Saturation  hysteresis  energy  density  (£^s)  7200  J/m3 

Saturation  applied  field  ( Hs )  75  kA/m 

Effective  relative  magnetic  permeability  (pr)  8.0 


H  ( kA/m ) 

Figure  12.  DC  Saturation  Hysteresis  Curve  for  Nickel  (from  Suwanwatana  [28]). 
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Frequency  (Hz) 

Figure  13.  Heating  Rate  as  a  Function  of  Frequency  for  Nickel  Particles  of  Various  Volume 
Fractions  in  a  Polymer  Film. 

From  Figure  7,  we  can  see  that  for  this  value  of  Ti  eddy  heating  becomes  important  for  parti¬ 
cles  in  the  1-10  /j,m  size  range  for  the  range  of  frequencies  we  are  considering.  Eddy  current 
heating  may  be  important  for  particles  which  are  larger  than  this  size.  Therefore,  to  ensure 
Curie  temperature  control,  the  particles  should  be  kept  smaller  than  1  fxm. 

Using  the  hysteresis  energy  density  approximation,  equation  35,  with  the  values  of  Ehs 
and  Hs  from  Table  2  and  a  peak  applied  field  magnitude  of  Hp  =  3500  A/m,  the  hysteresis 
power  density  for  nickel  is  15.7  J/m3.  Combining  this  value  with  equation  58  and  the  material 
properties  from  Table  2,  the  heating  rate  for  a  particle  /  polymer  film  is 

^  =  /(1.57  x  10-5)  (2.83  +  1/vf)-1 .  (65) 

Figure  13  shows  the  characteristic  heating  rate  as  a  function  of  frequency  for  various  volume 
fractions.  Figure  14  shows  the  same  information,  but  in  terms  of  volume  fractions  for  various 
heating  rates.  Figure  15  shows  heating  rates  as  a  function  of  particle  volume  fraction  for  various 
frequencies.  Note  that  the  frequency  must  be  at  least  1  MHz  to  get  reasonable  heating  rates. 
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Figure  14.  Nickel  Particle  Volume  Fraction  as  a  Function  of  Frequency  for  Various  Film 
Heating  Rates. 


Figure  15.  Heating  Rates  as  a  Function  of  Nickel  Particle  Volume  Fraction  for  Various  Fre¬ 
quencies. 


Table  3.  Strontium  Ferrite  Thermophysical  and  Magnetic  Properties.  All  values  above  hori¬ 
zontal  line  from  McCurrie  [31]  except  specific  heat,  which  is  based  on  typical  ferrite  values. 
All  values  below  line  determined  graphically  from  Figure  1 6. 


Property 

Value 

Density  (p) 

4800  kg/m3 

Specific  heat  (cp) 

750  J/kg  K 

Electrical  resistivity  (pa) 

104f2m 

Curie  temperature 

450°  C 

Saturation  hysteresis  energy  density  (Ehs) 

560  kJ/m3 

Saturation  applied  field  (Hs) 

600  kA/m 

Effective  relative  magnetic  permeability  (pr) 

1.7 

4.3  Strontium  Ferrite  (Hard  Ferrite).  Strontium  ferrite  (SrFe12019)  is  one  of  the 
most  common  commercial  permanent  magnet  materials.  It  is  cheap,  readily  available  as  a 
powder,  and  has  a  high  coercivity  typical  of  permanent  magnet  materials.  The  relevant  thermo¬ 
physical  and  magnetic  properties  for  a  typical  strontium  ferrite  are  shown  in  Table  3.  Figure  16 
shows  a  DC  saturation  hysteresis  curve  for  Triton  C  strontium  ferrite  [30],  developed  in  1996 
under  a  U.S.  Army  Small  Business  Innovative  Research  (SBIR)  program  with  Triton  Systems, 
Inc.  (Chelmsford,  MA)  for  the  Weapons  and  Materials  Research  Directorate,  U.S.  Army  Re¬ 
search  Laboratory.  The  magnetic  properties  in  Table  3  were  calculated  from  this  DC  hysteresis 
curve.  Ehs  is  the  area  of  the  loop,  Hs  is  applied  field  strength  necessary  to  mostly  saturate  the 
material,  and  p,r  is  measured  graphically,  as  suggested  in  section  2.3.  Note  that  the  field  nec¬ 
essary  to  reach  saturation  is  approximately  two  orders  of  magnitude  greater  than  the  fields  we 
expect  to  generate  in  our  induction  devices.  The  jj,r  value  of  1 .7  suggests  that  at  our  expected 
maximum  applied  field  strength,  3500  A/m,  the  material  magnetization  will  be  Bp  =  0.003  T. 

Using  equation  56  and  the  values  of  pa  and  j2r  from  Table  3,  the  penetration  thickness 
material  parameter  T2  has  a  value  for  strontium  ferrite  of 

T2  =  76.7  (66) 

From  Figure  8,  we  can  see  that  this  value  of  T2  does  not  even  appear  on  the  plot.  For  example, 
at  a  frequency  of  100  MHz,  shielding  effects  become  important  when  D0  =70  cm,  which  is 
much  larger  than  the  particles  we  will  be  using. 
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Figure  16.  DC  Saturation  Hysteresis  Curve  for  Strontium  Ferrite  (from  Stark  [30]). 

Using  equation  36  and  the  values  of  Ehs  and  Hs  from  Table  3,  we  calculate  ph  =  2.2. 
Substituting  this  value  with  the  values  of  pn  and  pr  into  equation  53,  the  heating  ratio  parameter 
T  i  has  a  value  for  strontium  ferrite  of 

Ti  =  2.4  x  10-4  f)_1m_1.  (67) 

This  value  of  Yi  does  not  even  appear  in  Figure  7.  For  example,  at  a  frequency  of  100  MHz, 
eddy  heating  becomes  important  when  D0  =  57  cm.  Therefore  Curie  temperature  control 
should  not  be  a  problem  for  strontium  ferrites. 

Using  the  hysteresis  energy  density  approximation,  equation  35,  with  the  values  of  Ehs 
and  Hs  from  Table  3  and  a  peak  applied  field  magnitude  of  Hp  =  3500  A/m,  the  hysteresis 
power  density  for  strontium  ferrite  is  19.1  J/m3.  Combining  this  value  with  equation  58  and 
the  material  properties  from  Table  3,  the  heating  rate  for  a  particle  /  polymer  film  is 

^  3  /(1.91  x  IQ'5)  (2.60  +  1/vj)-1 .  (68) 
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Frequency  (Hz) 

Figure  17.  Heating  Rate  as  a  Function  of  Frequency  for  Strontium  Ferrite  Particles  of  Various 
Volume  Fractions  in  a  Polymer  Film. 

Figure  17  shows  the  characteristic  heating  rate  as  a  function  of  frequency  for  various  volume 
fractions.  Figure  18  shows  the  same  information,  but  in  terms  of  volume  fractions  for  various 
heating  rates.  Figure  19  shows  heating  rates  as  a  function  of  particle  volume  fraction  for 
various  frequencies.  The  heating  rates  are  similar  to  those  observed  in  nickel,  with  frequencies 
of  1  MHz  or  higher  required  to  reach  l°C/s. 

4.4  Nickel  Zinc  Ferrite  (Soft  Ferrite).  Nickel  zinc  ferrite  (Nii_:rZnxFe204)  is  a 
soft  ferrite  which  is  used  extensively  in  transformer  cores,  mostly  due  to  its  high  permeability 
and  “low”  hysteresis  losses.  NiZn  ferrite  was  chosen  for  this  study  as  a  representative  soft 
ferrite  because  its  Curie  temperature  can  be  customized  from  100°-600°C  by  varying  the  ratio 
of  nickel  to  zinc  in  the  material  as  shown  in  Figure  20  [32].  Snelling  [33]  provides  extensive 
property  listings  and  hysteresis  curves  for  various  NiZn  ferrite  grades.  We  will  choose  a  grade 
of  NiZn  with  a  large  coercivity,  grade  B5  from  Snelling.  The  relevant  thermophysical  and 
magnetic  properties  for  a  typical  NiZn  ferrite  are  given  in  Table  4.  Figure  21  shows  a  DC 
saturation  hysteresis  loop  for  this  material.  The  magnetic  properties  in  Table  4  were  calculated 
from  this  DC  hysteresis  curve.  Ehs  is  the  area  of  the  loop,  Hs  is  applied  field  strength  necessary 
to  mostly  saturate  the  material,  and  jlr  is  measured  graphically,  as  suggested  in  section  2.3.  The 
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Figure  18.  Strontium  Ferrite  Particle  Volume  Fraction  as  a  Function  of  Frequency  for  Various 
Film  Heating  Rates. 


Figure  19.  Heating  Rates  as  a  Function  of  Strontium  Ferrite  Particle  Volume  Fraction  for 
Various  Frequencies. 


Table  4.  NiZn  Ferrite  Thermophysical  and  Magnetic  Properties.  All  values  above  horizontal 
line  from  Snelling  [33].  All  values  below  line  determined  graphically  from  Figure  21. 


Property 

Value 

Density  ( p ) 

4500  kg/m3 

Specific  heat  (cp) 

750  J/kg  K 

Electrical  resistivity  ( pa ) 

103fim 

Curie  temperature 

1 00-600°  C 

Saturation  hysteresis  energy  density  (Ehs) 

560  J/m3 

Saturation  applied  field  ( Hs ) 

4.0  kA/m 

Effective  relative  magnetic  permeability  (jlr) 

37.2 

fj,r  value  of  37.2  suggests  that  at  our  expected  maximum  applied  field  strength,  3500  A/m,  the 
material  magnetization  will  be  Bp  =  0.16  T. 

Using  equation  56  and  the  values  of  pa  and  jxr  from  Table  4,  the  penetration  thickness 
material  parameter  Y2  has  a  value  for  NiZn  ferrite  of 


T2  =  5.18  H1/2m1/2.  (69) 

From  Figure  8,  we  can  see  that  for  this  value  of  T2  shielding  effects  should  not  be  important. 
For  example,  at  100  MHz,  shielding  should  not  be  important  if  the  particle  diameters  are  less 
than  4  cm. 

Using  equation  36  and  the  values  of  Ehs  and  Hs  from  Table  4,  we  calculate  Jxh  -  28.9. 
Substituting  this  value  with  the  values  of  pa  and  Jxr  into  equation  53,  the  heating  ratio  parameter 
Ti  has  a  value  for  NiZn  ferrite  of 


Ti  =  5.0  x  10~2  ft-1m-1.  (70) 

This  value  of  Ti  does  not  even  appear  in  Figure  7.  For  example,  at  a  frequency  of  100  MHz, 
eddy  heating  becomes  important  when  D0  =4  cm.  Therefore  Curie  temperature  control  should 
not  be  a  problem  for  NiZn  ferrites. 

Using  the  hysteresis  energy  density  approximation,  equation  35,  with  the  values  of  Ehs  and 
Hs  from  Table  4  and  a  peak  applied  field  magnitude  of  Hp  =  3500  A/m,  the  hysteresis  power 
density  for  NiZn  ferrite  is  429  J/m3.  Combining  this  value  with  equation  58  and  the  material 
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Figure  20.  Curie  Temperature  as  a  Function  of  Composition  for  NiZn  Ferrite  (from  McCurrie 
[32]). 


Figure  21.  DC  Saturation  Hysteresis  Curve  for  NiZn  Ferrite  (from  Snelling  [33]). 
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Figure  22.  Heating  Rate  as  a  Function  of  Frequency  for  NiZn  Ferrite  Particles  of  Various 
Volume  Fractions  in  a  Polymer  Film. 

properties  from  Table  4,  the  heating  rate  for  a  particle  /  polymer  film  is 

dT 

—  “  /( 4.29  x  10"4)  (2.38  +  l/vj)-1 .  (71) 

Figure  22  shows  the  characteristic  heating  rate  as  a  function  of  frequency  for  various  volume 
fractions.  Figure  23  shows  the  same  information,  but  in  terms  of  volume  fractions  for  various 
heating  rates.  Figure  24  shows  heating  rates  as  a  function  of  particle  volume  fraction  for  various 
frequencies.  Note  that  reasonable  heating  rates  are  possible  at  frequencies  of  100  kHz.  These 
heating  rates  are  roughly  one  order  of  magnitude  greater  than  the  heating  rates  calculated  for 
nickel  and  strontium  ferrite.  At  10  MHz  and  vj  =  10%,  heating  rates  exceeding  100°  C/s  are 
possible. 

5.  Analysis 

5.1  Accuracy  of  Models.  Before  making  conclusions  and  recommendations  con¬ 
cerning  particle-based  induction  heating,  we  must  discuss  the  accuracy  of  the  calculations  per¬ 
formed  in  this  report.  The  eddy  current  model  (section  2.2),  given  accurate  material  properties. 
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Figure  23.  NiZn  Ferrite  Particle  Volume  Fraction  as  a  Function  of  Frequency  for  Various  Film 
Heating  Rates. 


Figure  24.  Heating  Rates  as  a  Function  of  NiZn  Ferrite  Particle  Volume  Fraction  for  Various 
Frequencies. 


should  be  very  accurate.  The  assumptions  of  uniform  field  and  properties  within  the  particle 
are  certainly  valid  for  the  small  particles  we  expect  to  be  using.  However,  because  we  do 
not  know  the  exact  material  permeability  and  resistivity  at  every  temperature,  frequency,  and 
applied  field  strength,  the  quantitative  accuracy  of  the  results  cannot  be  guaranteed.  We  can 
instead  expect  accuracy  within  one  or  two  orders  of  magnitude. 

The  hysteresis  heating  model  (section  2.3)  should  be  very  accurate,  given  accurate  hys¬ 
teresis  energy  density  functions.  For  the  case  of  NiZn  ferrite,  we  are  given  a  DC  hysteresis 
curve  that  is  within  our  applied  field  capabilities.  Therefore,  at  low  frequencies,  the  calculation 
should  be  extremely  accurate.  The  hysteresis  loop,  however,  may  change  dramatically  as  the 
frequency  is  increased,  so  a  high  level  of  accuracy  at  higher  frequencies  is  not  guaranteed.  For 
the  nickel  and  strontium  ferrite  cases,  our  hysteresis  heating  calculations  are  dependent  on  the 
estimation  minor  hysteresis  loops  based  on  a  saturation  hysteresis  measurement.  The  loop  self¬ 
similarity  approach  which  we  have  taken  is  the  a  reasonable  way  to  perform  these  calculations, 
but  the  accuracy  of  this  approach,  even  within  a  few  orders  of  magnitude,  is  not  guaranteed 
until  it  is  tested  experimentally. 

The  penetration  thickness  model  (section  2.4)  should  be  very  accurate,  although  we  are 
again  limited  by  the  accuracy  of  the  magnetic  permeability  values. 

The  induction  head  model  (section  2.5)  should  be  accurate  for  the  particular  geometry 
given.  However  it  must  be  emphasized  that  the  magnetic  field  strength  produced  by  a  piece  of 
equipment  is  determined  by  the  power  supply  specifications  and  the  geometry  of  the  induction 
coil.  We  have  performed  calculations  for  one  typical  coil  and  power  supply,  but  many  other 
options  are  available. 

The  heat  transfer  model  (section  2.6)  should  be  extremely  accurate,  and  the  thermophys¬ 
ical  material  properties  used  in  the  calculations  are  known  to  a  high  degree  of  accuracy.  The 
calculations  of  section  4.1  also  show  that,  as  long  as  the  particles  are  kept  to  a  1-10  [i m  size, 
the  uniform  heating  assumption  should  be  valid.  This  model  however  only  predicts  the  tem¬ 
perature  change  for  an  insulated  film.  In  joining  applications,  the  film  is  sandwiched  between 
two  other  non-heat  generating  composite  adherends.  The  adherends  act  as  heat  sinks,  drawing 
thermal  energy  away  from  the  bondline.  Therefore  the  actual  heating  rates  produced  in  joining 
experiments  will  be  less  than  those  predicted  by  this  model  and  measured  in  insulated  film 
experiments. 

In  summary,  the  fundamental  forms  of  all  of  the  models  are  well-grounded  in  the  physics 
of  heat  transfer  and  electromagnetism.  The  calculations  are  only  limited  by  the  accuracy  of 
the  material  properties,  especially  magnetic  material  properties.  However  we  should  expect 
that  all  calculations  are  accurate  to  within  one  or  two  orders  of  magnitude,  which  is  adequate 
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for  this  type  of  general  strategic  overview.  If  more  quantitative  modeling  is  to  be  done  in  the 
future,  significant  effort  must  be  invested  into  collecting  accurate  magnetic  material  data. 

5.2  Comparison  of  Materials.  The  calculations  of  section  4.2  show  that  nickel  par¬ 
ticles  will  provide  adequate  heating  rates  for  joining  of  composite  materials.  For  example, 
heating  rates  of  l°C/s  are  possible  at  a  frequency  of  1  MHz  and  a  volume  fraction  of  7%.  At 
10  MHz,  a  dispersion  with  7%  volume  fraction  would  heat  at  10° C/s.  However,  size  issues  may 
become  important.  For  example,  our  calculations  showed  that  at  high  frequencies  the  particles 
may  need  to  be  less  than  1  /j.m  in  order  to  prevent  shielding  and  eddy  current  effects.  If  our 
material  property  data  is  inaccurate,  these  maximum  particle  sizes  could  be  even  smaller.  As 
particle  sizes  approach  the  nanometer  size  range,  special  problems  arise.  These  difficulties  are 
discussed  in  section  5.3. 

The  calculations  of  section  4.3  show  that  strontium  ferrite  particles  will  provide  adequate 
heating  rates  for  joining  of  composite  materials  under  similar  conditions  as  nickel  particles. 
However,  unlike  nickel,  shielding  and  eddy  heating  effects  are  insignificant  for  any  particles 
less  than  1  cm  in  diameter.  This  difference  is  due  to  the  nine  orders  of  magnitude  difference 
between  the  resistivity  of  nickel  and  strontium  ferrite.  The  low  resistivity  of  ferrites  is  therefore 
extremely  beneficial.  Even  at  very  high  frequencies,  shielding  and  eddy  heating  effects  will  not 
constrain  the  maximum  particle  size,  and  we  can  avoid  any  potential  difficulties  associated  with 
nanometer-sized  particles. 

The  calculations  of  section  4.4  show  that  NiZn  ferrite  particles  will  provide  excellent  heat¬ 
ing  rates  for  joining  of  composite  materials,  even  at  moderate  frequencies  and  low  volume 
fractions.  For  example,  heating  rates  of  50°C/s  are  possible  at  a  frequency  of  1  MHz  and  a 
volume  fraction  of  15%.  NiZn  ferrites  also  have  the  high  resistivity  common  to  all  ferrites,  so 
that  shielding  and  eddy  heating  effects  are  never  important. 

A  comparison  of  the  important  behaviors  of  the  nickel,  strontium  ferrite,  and  NiZn  ferrite 
materials  is  given  in  Table  5.  From  the  table  it  is  clear  that  shielding  and  eddy  heating  effects 
may  limit  the  maximum  frequency  for  processing  using  nickel  particles,  but  no  practical  fre¬ 
quency  limit  exists  for  processing  using  ferrite  particles.  The  comparison  of  hysteresis  energy 
density  Eh  also  shows  why  NiZn  ferrite  provides  a  heating  rate  one  order  of  magnitude  greater 
than  nickel  and  strontium  ferrite,  since  heating  rate  is  proportional  to  Ek  (equation  58). 

This  hysteresis  energy  density  comparison  is  shown  graphically  in  Figure  25  for  an  alter¬ 
nating  field  of  maximum  strength  3500  A/m.  Because  the  applied  field  strength  is  similar  to 
the  Hs  of  the  NiZn  ferrite,  its  hysteresis  loop  is  utilized  fully.  Nickel  has  a  much  higher  Hs  and 
lower  p,r,  so  it  is  not  fully  magnetized  at  this  field  strength.  It  is  also  evident  that  the  energy 
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Table  5.  Performance  Comparison  of  Magnetic  Materials. 


Property 

Nickel 

Strontium  ferrite 

NiZn  ferrite 

Frequency  at  which  D0= 1  /im  particle 
has  significant  shielding 

68  MHz 

4.7  x  1010  GHz 

2.1  x  108  GHz 

Frequency  at  which  D0  —  1/j.m  particle 
has  significant  eddy  heating 

8.5  MHz 

3.3  x  1010  GHz 

1.6  x  10s  GHz 

Eh  for  Hp= 3500  A/m 

15.7  J/m3 

19.1  J/m3 

429  J/m3 

Figure  25.  Comparison  of  Hysteresis  Curves  of  Nickel,  NiZn  Ferrite,  and  Sr  Ferrite  in  an  AC 
Field  with  Maximum  Field  of  3500  A/m. 


39 


density  of  nickel  is  reduced  because  its  hysteresis  loop  is  not  very  “square”  (Figure  12).  Com¬ 
paring  the  hysteresis  curve  for  strontium  ferrite  in  Figure  25  with  its  saturation  hysteresis  curve 
(Figure  1 6)  dramatically  illustrates  the  ineffectiveness  of  using  hard  ferrites  (permanent  magnet 
materials)  for  heating  in  low-strength  fields.  Hard  ferrites  have  larger  DC  saturation  hysteresis 
loop  areas  than  soft  ferrites.  However,  the  coercivities  necessary  to  reach  saturation  in  these 
materials  (~  100  —  1000  kA/m)  is  many  orders  of  magnitude  greater  than  we  can  ever  apply 
with  practical  portable  equipment,  and  the  resulting  saturation  magnetization  is  not  necessarily 
any  greater  than  the  saturation  magnetization  of  soft  ferrites.  This  statement  is  equivalent  to 
the  observation  that  hard  ferrites  have  extremely  low  permeabilities.  Permanent  magnet  mate¬ 
rials  are  designed  to  resist  changes  in  magnetization  by  an  applied  magnetic  field.  Therefore, 
in  the  presence  of  modest  alternating  magnetic  fields,  the  resulting  change  in  magnetization  is 
extremely  low  and  hysteresis  heating  is  negligible. 

The  generalization  that  soft  ferrites  are  superior  to  hard  ferrites  is  not  entirely  true,  however. 
The  typical  coercivity  range  of  soft  ferrites  is  only  ~  1  —  100  A/m,  which  is  much  smaller  than 
the  applied  fields  our  equipment  should  be  capable  of.  We  will  be  able  to  easily  saturate  the 
material,  but  the  total  energy  required  to  do  so  is  not  significant.  The  most  efficient  heating 
instead  occurs  if  the  susceptor  material  has  a  coercivity  roughly  equal  to  our  maximum  applied 
field,  ~  1000  A/m,  as  was  the  case  for  the  particular  grade  of  NiZn  ferrite  which  was  used  for 
our  calculations.  However,  this  value  of  coercivity  is  unusually  high  for  a  soft  ferrite.  The  types 
of  materials  which  are  best  suited  for  our  heating  applications  would  probably  be  called  “semi- 
hard”  magnetic  materials.  Because  of  their  moderate  coercivities,  they  would  prove  inferior  to 
other  materials  for  most  conventional  hard  and  soft  magnetic  material  applications.  For  this 
reason  there  is  very  little  research  on  or  commercialization  of  these  materials. 

Although  the  coercivity  of  the  particular  NiZn  ferrite  grade  used  in  this  report  is  nearly 
ideal,  the  saturation  magnetization  is  actually  relatively  low,  around  0.15  T.  As  a  comparison, 
the  saturation  magnetization  of  elemental  iron  is  over  2  T.  Therefore  there  is  a  potential  order 
of  magnitude  gain  in  the  hysteresis  energy  density  if  the  proper  material  can  be  found. 

5.3  Size  Effects.  Using  very  small  particles  has  a  number  of  advantages.  For  a  par¬ 
ticulate  film  with  a  given  volume  fraction  of  particles,  the  uniformity  of  heat  generation  will 
be  higher  if  the  particles  are  smaller,  since  the  average  spacing  between  particles  is  reduced. 
Shielding  effects  and  eddy  current  heating,  phenomena  which  limit  the  applicability  of  our 
approach,  are  reduced  as  particle  size  is  decreased.  Smaller  particles  should  provide  less  of 
a  detriment  to  the  mechanical  properties  of  the  film  than  large  particles.  This  feature  is  espe¬ 
cially  important  for  composite  materials,  whose  micrometer-sized  microstructure  is  carefully 
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designed  for  optimum  macroscopic  properties.  Smaller  particles  also  should  have  less  of  an 
effect  on  the  effective  viscosity  of  the  film  during  processing,  as  compared  with  larger  particles. 

However,  as  particle  sizes  approach  the  nanometer  size  range,  special  problems  arise.  In 
order  to  produce  the  films,  the  particles  must  be  dispersed  in  a  polymer  matrix.  Melt  process¬ 
ing  is  the  most  likely  approach.  Nanometer-sized  particles  have  significant  inter-particle  forces 
which  may  make  dispersion  difficult.  For  metal  particles,  oxidation  is  an  additional  problem. 
The  extremely  high  surface  area  to  mass  ratio  of  nanometer-sized  particles  causes  rapid,  high- 
conversion  percentage  oxidation.  For  nickel,  the  resulting  oxide  is  not  magnetic  or  conductive, 
so  that  the  heating  ability  of  the  particles  is  greatly  reduced  or  even  eliminated.  Special  en¬ 
vironmental  handling  or  encapsulation  could  be  used  to  prevent  oxidation,  but  would  require 
additional  research  to  implement. 

Production  of  nanometer-sized  particles  is  also  difficult,  particularly  with  metallic  materi¬ 
als.  Micrometer-sized  particles  can  be  created  relatively  easily  using  mechanical  methods.  In 
contrast,  nanometer-sized  particles  require  either  extremely  precise  (and  expensive)  mechani¬ 
cal  equipment  or  special  techniques  such  as  sputtering  or  laser  ablation.  These  smaller  particle 
methods  also  typically  have  much  lower  yields  than  more  conventional  mechanical  techniques. 
Furthermore,  alloying  materials  is  difficult  at  such  small  size  scales,  because  the  individual 
constituents  tend  to  separate. 

6.  Conclusions 

6.1  Ideal  Material  and  Equipment  Characteristics.  From  our  calculations 
and  analysis,  we  can  now  state  the  ideal  material  and  equipment  characteristics  for  optimum 
induction  heating  performance  and  applicability.  We  seek  a  material  with  very  high  resistivity, 
in  order  to  minimize  shielding  effects  and  eddy  current  heating.  The  hysteresis  loop  in  the 
presence  of  a  high  frequency  (>  10  kHz)  and  moderate  amplitude  (~  3500  A/m)  field  should 
have  maximum  area.  In  terms  of  magnetic  material  properties,  high  hysteresis  losses  under 
such  conditions  requires  high  relative  magnetic  permeability  (>  100)  and  coercivities  in  the  1- 
10  kA/m  range.  Loop  “squareness”  is  also  beneficial.  We  are  ultimately  interested  in  polymer 
processing  between  100°  and  400°  C,  and  would  like  to  be  able  to  customize  our  processing 
temperature  within  ~  10°C.  Rather  than  try  to  find  25  different  materials,  each  with  a  Curie 
temperature  at  a  different  temperature,  we  would  rather  find  one  material  or  class  of  materials 
whose  Curie  temperature  can  be  precisely  controlled  through  alloying  or  processing  history. 
The  NiZn  ferrite  is  an  excellent  example  of  this,  with  its  Curie  temperature  very  precisely 
determined  by  the  relative  proportions  of  nickel  and  zinc  ferrite  in  the  structure. 


41 


The  size  of  the  particles  should  be  kept  small,  to  avoid  shielding  and  eddy  current  effects, 
to  minimize  particle  influence  on  bulk  mechanical  and  viscous  properties,  and  to  make  heating 
in  the  film  uniform.  However,  the  particles  should  be  kept  above  the  nanometer  size  range  to 
avoid  oxidation,  dispersion,  and  production  difficulties.  The  1-10  //m  size  range  seems  to  offer 
the  best  compromise  of  these  properties. 

Higher  frequency  induction  heating  equipment  provides  higher  heating  rates  in  susceptor 
materials,  until  shielding  effects  begin  to  limit  the  field  strength  in  the  particle.  For  1  fim 
diameter  conductive  particles,  this  upper  frequency  is  around  10  MHz.  For  the  insulating 
ferrites,  the  upper  frequency  is  much  higher.  It  is  also  possible  that  the  hysteresis  mechanisms 
in  the  particle  change  at  very  high  frequencies.  However,  until  such  limiting  mechanisms  are 
identified,  very  high  frequency  field  sources,  including  microwave  sources,  should  be  pursued. 

6.2  Summary.  Some  of  the  most  important  conclusions  from  this  study  are: 

•  Metallic  particles  produce  reasonable  heating  rates  but  are  susceptible  to  shielding  and 
eddy  heating  effects  at  higher  frequencies. 

•  Low  resistivity  materials,  such  as  ferrites,  are  not  susceptible  to  shielding  and  eddy  heat¬ 
ing  effects  at  practical  frequencies. 

•  Hard  magnetic  materials  have  generally  too  low  a  permeability  and  too  high  a  coercivity 
to  be  useful  for  hysteresis  heating. 

•  Soft  magnetic  materials  generally  have  adequate  permeability  to  be  useful  for  hysteresis 
heating,  but  often  have  too  low  of  a  coercivity. 

•  The  ideal  material  for  hysteresis  heating  applications  has  a  coercivity  in  a  range  which 
is  not  useful  for  most  other  magnetic  applications.  For  this  reason,  such  “semi-hard” 
materials  have  not  been  investigated  significantly. 

•  With  NiZn  ferrite  particles,  excellent  heating  rates  can  be  attained  with  available  equip¬ 
ment  (1-10  kW,  100  kHz-10  MHz)  and  low  particle  volume  fractions  (~  10%). 

•  Keeping  particle  sizes  in  the  1-10  fi m  range  should  provide  uniform  film  heating,  elimi¬ 
nate  shielding  and  eddy  effects,  and  allow  easy  particle  production. 

•  Accurate  mathematical  models  of  the  relevant  process  physics  for  the  particle  film  in¬ 
duction  heating  process  can  be  written.  However,  the  accuracy  of  their  implementation 
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is  seriously  limited  by  the  unavailability  of  the  proper  magnetic  property  data  in  the 
literature. 


7.  Recommendations  for  Future  Work 

Based  on  the  conclusions  of  this  report,  we  can  recommend  a  number  of  research  thrusts 
which  will  allow  for  the  implementation  and  sophistication  of  the  particle  induction  heating 
process.  These  research  areas,  which  are  discussed  in  the  following  sections,  include: 

•  Preliminary  induction  heating  experiments  with  NiZn  ferrite  particle  /  polymer  films  and 
adhesive  pastes, 

•  A  comprehensive  literature  and  commercial  search  to  find  candidate  materials  with  high 
potential  for  hysteresis  heating  applications, 

•  Design  and  implementation  of  a  device  for  direct  measurement  of  Eh, 

•  Experimentation  to  determine  the  feasibility  of  microwave  induction  processing, 

•  Literature  search  on  magnetic  resonance  phenomenon,  to  determine  the  feasibility  of 
improved  heating  through  precise  selection  of  source  frequency,  and 

•  Determination  of  the  potential  enhancement  of  hysteresis  energy  density  through  the  use 
of  nanometer-sized  particles. 

7.1  Preliminary  Experiments  with  NiZn  Ferrite.  The  calculations  of  section  4.4 
show  that  NiZn  ferrite  should  provide  excellent  heating  rates  in  magnetic  fields  of  100  kHz- 
1  MHz,  which  is  within  the  frequency  range  of  existing  equipment  at  the  U.S.  Army  Research 
Laboratory  and  the  University  of  Delaware  Center  for  Composite  Materials  (CCM).  NiZn  fer¬ 
rite  is  also  available  commercially,  is  inexpensive,  and  can  be  found  in  many  grades  (e.  g.  Stew¬ 
ard  USA,  Chatanooga,  TN;  Ceramic  Magnets,  Fairfield,  NJ;  Fair-Rite  Products,  Wallkill,  NY; 
Powder  Tech,  Valparaiso,  IN).  NiZn  ferrite  powder  should  be  acquired  and  simple  film  heating 
tests  performed.  Consolidating  a  film  /  powder  /  film  sandwich  in  a  hot  press  can  be  used  to 
manufacture  the  film.  The  film  can  be  insulated  and  placed  in  a  known  magnetic  field,  and  the 
temperature  measured  as  a  function  of  time. 
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7.2  Magnetic  Materials  Literature  Search.  Extensive  literature  exists  on  mag¬ 
netic  materials  and  their  properties.  A  literature  search  needs  to  be  performed  to  find  general 
classes  of  materials  and  compositional  families  which  have  the  optimal  properties  described  in 
section  6.1.  This  search  should  include,  but  not  be  limited  to,  scientific  journals,  patent  dis¬ 
closures,  and  commercial  databases.  It  will  also  be  important  to  determine  general  effects  of 
specific  elemental  additives  or  processing  conditions  on  relevant  magnetic  properties,  which 
permit  exploration  of  previously  untested  materials. 

7.3  Direct  Measurement  of  Eh •  In  order  to  assess  and  optimize  the  heating  po¬ 
tential  of  a  material,  its  hysteresis  energy  density  Eh  needs  to  be  measured  as  a  function  of 
frequency,  applied  field  amplitude,  and  temperature.  Equipment  exists  for  measuring  initial 
permeability  and  loss  tangents  as  a  function  of  frequency,  and  DC  hysteresis  loops  as  a  func¬ 
tion  of  applied  field  amplitude  and  temperature.  However,  all  of  these  measurements  only 
allow  us  to  partially  infer  Eh  for  a  given  set  of  conditions. 

A  tool  is  needed  which  can  directly  measure  £/,asa  function  of  frequency,  applied  field 
amplitude,  and  temperature.  The  applied  field  strength  must  be  variable  up  to  10  kA/m,  and 
have  a  frequency  range  of  10  kHz- 100  MHz.  The  measurements  should  only  require  a  small 
sample  mass,  so  that  expensive  or  difficult  to  attain  materials  can  be  tested.  The  device  should 
also  have  temperature  measurement  and  control  capabilities. 

An  extensive  search  needs  to  be  undertaken  to  determine  if  such  a  device  already  exists. 
Devices  which  measure  loss  tangents  are  candidates,  as  well  as  transformer  core  loss  measure¬ 
ment  equipment. 

If  existing  devices  cannot  be  found  that  meet  our  requirements,  new  equipment  needs  to 
be  designed  and  constructed.  The  device  needs  to  generate  an  oscillating  magnetic  field  in  a 
small  sample  of  material.  The  power  consumed  can  be  measured  electrically,  in  terms  of  power 
supply  power,  or  thermally,  for  example  through  measurement  of  temperature  rise  of  a  nitrogen 
stream  flowing  past  the  sample.  The  unit  must  also  have  an  integrated  temperature  control  and 
inert  atmosphere  capability.  Note  that  in  both  of  these  scenarios,  no  magnetic  measurements 
need  to  made.  We  are  only  measuring  power  consumption  in  a  mass  of  material.  Construction 
of  such  a  device  should  be  possible  with  reasonable  resources. 

7.4  Microwave  Processing.  According  to  our  models,  the  heating  power  in  the  film 
increases  proportionally  with  frequency.  Our  calculations  also  show  that  with  ferrites  shield¬ 
ing  does  not  become  important  until  very  high  frequencies.  Therefore  microwave  frequency 
magnetic  waves  (/  >  1  GHz)  should  provide  extremely  efficient  heating.  It  is  possible  that  at 
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such  high  frequencies,  the  hysteresis  heat  generation  mechanisms  change.  Experiments  will 
provide  qualitative  evidence  of  the  feasibility  of  microwave  heating.  Candidate  equipment 
includes  standard  1  kW  2.45  GHz  microwave  ovens,  or  the  6  kW  2.45  GHz  multi-mode  mi¬ 
crowave  chamber  at  CCM.  Based  on  our  earlier  calculations,  NiZn  ferrite  particles  could  be 
used  as  the  susceptor  material. 

7.5  Magnetic  Resonance.  All  of  the  calculations  in  this  report  assume  that  the  hys¬ 
teresis  energy  density  of  the  susceptor  materials  does  not  change  with  frequency.  However, 
at  very  high  frequencies,  this  assumption  is  not  correct.  The  material  does  not  have  sufficient 
time  to  respond  to  the  applied  magnetic  field,  and  so  it  does  not  fully  magnetize.  We  expect 
that  this  effect  will  lower  the  energy  density  of  the  material  as  frequency  increases.  However, 
at  these  high  frequencies,  low-amplitude  magnetic  losses  become  significant.  At  a  critical  fre¬ 
quency  the  material  is  in  a  state  of  magnetic  resonance  [34],  where  low -amplitude  losses  in 
the  material  are  maximized.  Magnetic  resonance  frequencies  are  typically  in  the  100  MHz- 
10  GHz  range.  The  existence  of  this  resonant  behavior  implies  that  it  may  be  advantageous 
to  “tune”  the  applied  field  frequency  to  a  specific  frequency,  rather  than  trying  to  maximize 
power  through  maximization  of  source  frequency.  It  is  not  clear,  however,  whether  magnetic 
resonance  can  produce  energy  densities  comparable  to  lower  frequency,  high  amplitude  hys¬ 
teresis  mechanisms.  A  careful  investigation,  primarily  through  literature  investigation,  needs 
to  be  performed  to  determine  the  importance  of  magnetic  resonance  effects. 

7.6  Submicron  Particle  Issues.  The  magnetic  properties  of  materials  become  size- 
dependent  as  the  particle  size  approaches  the  nanometer  size  scale.  It  is  possible  that  reducing 
particle  size  to  the  nanometer  size  scale  could  result  in  significant  enhancement  of  the  hys¬ 
teresis  energy  density  of  a  material.  Eh  needs  to  be  measured  as  a  function  of  particle  size  to 
determine  if  this  effect  is  significant. 

Many  other  issues  associated  with  nanoparticles  exist,  including  particle  production,  dis¬ 
persion,  and  oxidation  prevention.  However,  all  of  these  potential  problems  can  be  avoided 
by  using  micron-sized  particles.  In  addition,  this  report  finds  no  overwhelming  advantages 
to  using  nanometer-sized  particles.  Therefore,  the  investigation  of  the  production,  dispersion, 
and  oxidation  of  nanometer  particles  should  not  be  pursued  until  clear  benefits  of  nanometer 
particles  are  identified. 
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4  DIRECTOR 

US  ARMY  CECOM 
NIGHT  VISION  & 

ELECTRONIC  SENSORS  DIR 

AMSEL  RD  NV  CM  CCD 

R  ADAMS 

R  MCLEAN 

A  YINGST 

AMSEL  RD  NV  VISP 

E  JACOBS 

10221  BURBECK  RD 

FT  BELVOIR  VA  22060-5806 

2  US  ARMY  CORPS  OF  ENGINEERS 
CERDC 
T  LIU 
CEW  ET 
TTAN 

20  MASS  AVE  NW 
WASHINGTON  DC  20314 


HQ  IOC  TANK 
AMMUNITION  TEAM 
AMSIO  SMT 
R  CRAWFORD 
W  HARRIS 

ROCK  ISLAND  IL  61299-6000 

DIRECTOR 
US  ARMY  AMCOM 
SFAE  AV  RAM  TV 
D  CALDWELL 
BLDG  5300 

REDSTONE  ARSENAL  AL 
35898 

COMMANDER 
US  ARMY  AMCOM 
AVIATION  APPLIED  TECH  DIR 
J  SCHUCK 

FT  EUSTIS  VA  23604-5577 

US  ARMY  CERL 
RLAMPO 

2902  NEWMARK  DR 
CHAMPAIGN  IL  61822 


1  US  ARMY  COLD  REGIONS 
RSCH  &  ENGRNG  LAB 
P  DUTTA 
72  LYME  RD 
HANOVER  NH  03755 

1  SYSTEM  MANAGER  ABRAMS 

ATZK  TS 
LTC  JHNUNN 
BLDG  1002  RM  110 
FT  KNOX  KY  40121 

1  CHIEF 

USAIC 
ATZB  COM 
LTC  T  J  CUMMINGS 
FT  BENNING  GA  3 1 905-5800 

1  USA  SBCCOM  PM  SOLDIER  SPT 

AMSSB  PM  RSS  A 
J  CONNORS 
KANSAS  ST 

NATICK  MA  01760-5057 
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8  DIRECTOR 

US  ARMY  NATIONAL  GROUND 
INTELLIGENCE  CTR 
D  LEITER 
M  HOLTUS 
M  WOLFE 
S  MINGLEDORF 
J  GASTON 

W  GST  ATTENB  AUER 
R  WARNER 
J  CRIDER 

220  SEVENTH  ST  NE 
CHARLOTTESVILLE  VA  22091 

6  US  ARMY  SBCCOM 

SOLDIER  SYSTEMS  CENTER 
BALLISTICS  TEAM 
J  WARD 

MARINE  CORPS  TEAM 

J  MACKIEWICZ 

BUS  AREA  ADVOCACY  TEAM 

W  HASKELL 

SSCNCWST 

W  NYKVIST 

T  MERRILL 

S  BEAUDOIN 

KANSAS  ST 

NATICK  MA  01760-5019 

9  US  ARMY  RESEARCH  OFC 
ACROWSON 

J CHANDRA 
H EVERETT 
J PRATER 
R  SINGLETON 
G ANDERSON 
D  STEPP 
D  KISEROW 
J CHANG 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

1  NAVAL  SURFACE  WARFARE  CTR 
DAHLGREN  DIV  CODE  G06 
DAHLGREN  VA  22448 


NO.  OF 

COPIES  ORGANIZATION 

1  NAVAL  SURFACE  WARFARE  CTR 

TECH  LIBRARY  CODE  323 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

3  NAVAL  RESEARCH  LAB 

I WOLOCK  CODE  6383 
R  BADALIANCE  CODE  6304 
L  GAUSE 

WASHINGTON  DC  20375 

1  NAVAL  SURFACE  WARFARE  CTR 
CRANE  DIVISION 

M  JOHNSON  CODE  20H4 
LOUISVILLE  KY  40214-5245 

2  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
CARDEROCK  DIVISION 
R  PETERSON  CODE  2020 
M  CRITCHFIELD  CODE  1730 
BETHESDA  MD  20084 

2  NAVAL  SURFACE  WARFARE  CTR 

U  SORATHIA 
C  WILLIAMS  CD  6551 
9500  MACARTHUR  BLVD 
WEST  BETHESDA  MD  20817 

1  DAVID  TAYLOR  RESEARCH  CTR 
SHIP  STRUCTURES 

&  PROTECTION 
DEPT  CODE  1702 
BETHESDA  MD  20084 

2  DAVID  TAYLOR  RESEARCH  CTR 
R  ROCKWELL 

W  PHYILLAIER 
BETHESDA  MD  20054-5000 

1  OFC  OF  NAVAL  RESEARCH 

D  SIEGEL  CODE  351 
800  N  QUINCY  ST 
ARLINGTON  VA  22217-5660 
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8  NAVAL  SURFACE  WARFARE  CTR 
J  FRANCIS  CODE  G30 
D  WILSON  CODE  G32 
R  D  COOPER  CODE  G32 
J  FRAYSSE  CODE  G33 
E  ROWE  CODE  G33 
T  DURAN  CODE  G33 
L  DE  SIMONE  CODE  G33 
R  HUBBARD  CODE  G33 
DAHLGREN  VA  22448 

1  NAVAL  SEA  SYSTEMS  CMD 

D  LIESE 

2531  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22242-5160 

1  NAVAL  SURFACE  WARFARE  CTR 

M  LACY  CODE  B02 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

1  OFC  OF  NAVAL  RES 
J  KELLY 

800  NORTH  QUINCEY  ST 
ARLINGTON  VA  22217-5000 

2  NAVAL  SURFACE  WARFARE  CTR 
CARDEROCK  DIVISION 

R  CRANE  CODE  2802 
C  WILLIAMS  CODE  6553 
3A  LEGGETT  CIR 
BETHESDA  MD  20054-5000 

1  EXPEDITIONARY  WARFARE  DIV  N85 

FSHOUP 

2000  NAVY  PENTAGON 
WASHINGTON  DC  20350-2000 

1  AFRL  MLBC 

2941  PSTRM  136 

WRIGHT  PATTERSON  AFB  OH 

45433-7750 

1  AFRL  MLSS 

R  THOMSON 
2179  12TH  STRM  122 
WRIGHT  PATTERSON  AFB  OH 
45433-7718 


NO.  OF 

COPIES  ORGANIZATION 

2  AFRL 

F ABRAMS 

J  BROWN 

BLDG  653 

2977  P  ST  STE  6 

WRIGHT  PATTERSON  AFB  OH 

45433-7739 

1  AFRL  MLS  OL 

L  COULTER 
7278  4TH  ST 
BLDG  100  BAY  D 
HILL  AFB  UT  84056-5205 

1  OSD 

JOINT  CCD  TEST  FORCE 
OSD  JCCD 
R  WILLIAMS 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  29180-6199 

1  DEFENSE  NUCLEAR  AGENCY 

INNOVATIVE  CONCEPTS  DIV 
6801  TELEGRAPH  RD 
ALEXANDRIA  VA  22310-3398 

1  WATERWAYS  EXPERIMENT 
D SCOTT 

3909  HALLS  FERRY  RD  SC  C 
VICKSBURG  MS  39180 

3  DARPA 

M  VANFOSSEN 
S  WAX 

L  CHRISTODOULOU 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

2  FAA 

TECH  CENTER 
D  OPLINGER  AAR  431 
P  SHYPRYKEVICH  AAR  431 
ATLANTIC  CITY  NJ  08405 

2  SERDP  PROGRAM  OFC 

PMP2 

C  PELLERIN 
B  SMITH 

901  N  STUART  ST  STE  303 
ARLINGTON  VA  22203 
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NO.  OF  NO.  OF 

COPIES  ORGANIZATION  COPIES 


1  FAA  7 

MIL  HDBK  17  CHAIR 
LILCEWICZ 
1601  LIND  AVESW 
ANM  115N 
RESTON  VA  98055 

1  US  DEPT  OF  ENERGY 

OFC  OF  ENVIRONMENTAL 
MANAGEMENT 
P  RITZCOVAN 

19901  GERMANTOWN  RD  1 

GERMANTOWN  MD  20874-1928 

1  DIRECTOR 

LLNL 

F  ADDESSIO  MSB216 

PO  BOX  1633  3 

LOS  ALAMOS  NM  87545 

5  DIRECTOR 

LLNL 

R  CHRISTENSEN 
S  DETERESA 
F  MAGNESS 
M  FINGER  MS  313 
M  MURPHY  L  282 

PO  BOX  808  3 

LIVERMORE  CA  94550 

1  OAK  RIDGE  NATIONAL 

LABORATORY 
RM  DAVIS 
PO  BOX  2008 

OAK  RIDGE  TN  37831-6195  1 

1  OAK  RIDGE  NATIONAL 

LABORATORY 

C  EBERLE  MS  8048  1 

PO  BOX  2009 

OAK  RIDGE  TN  37831 

1  OAK  RIDGE  NATIONAL  1 

LABORATORY 
CD  WARREN  MS  8039 
PO  BOX  2009 
OAK  RIDGE  TN  37922 


ORGANIZATION 

NIST 

RPARNAS 

JDUNKERS 

M  VANLANDINGHAM  MS  8621 
J  CHIN  MS  8621 
D  HUNSTON  MS  8543 
J  MARTIN  MS  8621 
D  DUTHINH  MS  8611 
100  BUREAU  DR 
GAITHERSBURG  MD  20899 

HYDROGEOLOGIC  INC 
SERDP  ESTCP  SPT  OFC 
S  WALSH 

1155  HERNDON  PKWY  STE  900 
HERNDON  VA  20170 

DIRECTOR 

SANDIA  NATIONAL  LABS 
APPLIED  MECHANICS  DEPT 
DIV  8241 
I HANDROCK 
Y  R  KAN 
J  LAUFFER 
PO  BOX  969 

LIVERMORE  CA  94550-0096 

NASA  LANGLEY  RSCH  CTR 
AMSRLVS 
WELBERMS  266 
F  BARTLETT  JR  MS  266 
G  FARLEY  MS  266 
HAMPTON  VA  23681-0001 

NASA  LANGLEY  RSCH  CTR 
T  GATES  MS  188E 
HAMPTON  VA  23661-3400 

USDOT  FEDERAL  RAILRD 
M  FATEH  RDV  31 
WASHINGTON  DC  20590 

FHWA 

EMUNLEY 

6300  GEORGETOWN  PIKE 
MCLEAN  VA  22101 
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ORGANIZATION 


NO.  OF 
COPIES 


1  CENTRAL  INTLLGNC  AGNCY 

OTI  WDAG  GT 
W  L  WALTMAN 
PO  BOX  1925 
WASHINGTON  DC  20505 

1  MARINE  CORPS  INTLLGNC  ACTVTY 

D  KOSITZKE 

3300  RUSSELL  RD  STE  250 
QUANTICO  VA  22134-501 1 

1  DIRECTOR 

NATIONAL  GRND  INTLLGNC  CTR 
IANG  TMT 
220  SEVENTH  ST  NE 
CHARLOTTESVILLE  VA 
22902-5396 

1  DIRECTOR 

DEFENSE  INTLLGNC  AGNCY 
TA  5 

K  CRELLING 
WASHINGTON  DC  203 1 0 

1  GRAPHITE  MASTERS  INC 

J  WILLIS 

3815  MEDFORD  ST 

LOS  ANGELES  CA  90063-1900 

1  ADVANCED  GLASS  FIBER  YARNS 

T  COLLINS 

281  SPRING  RUN  LANE  STE  A 
DOWNINGTON  PA  19335 

1  COMPOSITE  MATERIALS  INC 

D  SHORTT 
19105  63  AVENE 
PO  BOX  25 

ARLINGTON  WA  98223 

1  COMPOSITE  MATERIALS  INC 

R  HOLLAND 
11  JEWEL  CT 
ORINDA  CA  94563 

1  COMPOSITE  MATERIALS  INC 

C  RILEY 

14530  S  ANSON  AVE 
SANTA  FE  SPRINGS  CA  90670 


2  COMPOSIX 

D BLAKE 
L  DIXON 
120  O  NEILL  DR 
HEBRUN  OHIO  43025 

4  CYTEC  FIBERITE 

R  DUNNE 
D  KOHLI 
M  GILLIO 
RMAYHEW 
1300  REVOLUTION  ST 
HAVRE  DE  GRACE  MD  21078 

2  SIMULA 

J  COLTMAN 
RHUYETT 
10016  S51ST  ST 
PHOENIX  AZ  85044 

1  SIOUX  MFG 
BKRIEL 
PO  BOX  400 

FT  TOTTEN  ND  58335 

2  PROTECTION  MATERIALS  INC 
M  MILLER 

F  CRILLEY 
14000  NW  58  CT 
MIAMI  LAKES  FL  33014 

3  FOSTER  MILLER 
J  J  GASSNER 

M  ROYLANCE 
WZUKAS 
195  BEAR  HILL  RD 
WALTHAM  MA  02354-1 196 

1  ROM  DEVELOPMENT  CORP 
ROMEARA 

136  SWINEBURNE  ROW 
BRICK  MARKET  PLACE 
NEWPORT  RI  02840 

2  TEXTRON  SYSTEMS 
T  FOLTZ 

M  TREASURE 
201  LOWELL  ST 
WILMINGTON  MA  08870-2941 
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1  JPS  GLASS 

L CARTER 
PO  BOX  260 
SLATER  RD 
SLATER  SC  29683 

1  O  GARA  HESS  &  EISENHARDT 
M  GILLESPIE 

9113  LESAINTDR 
FAIRFIELD  OH  45014 

2  MILLIKEN  RSCH  CORP 
H  KUHN 

M  MACLEOD 
PO  BOX  1926 
SPARTANBURG  SC  29303 

1  CONNEAUGHT  INDUSTRIES  INC 
J  SANTOS 

PO  BOX  1425 
COVENTRY  RI  02816 

2  BATTELLE  NATICK  OPNS 
J  CONNORS 
BHALPIN 

209  W  CENTRAL  ST  STE  302 
NATICK  MA  01760 

1  ARMTEC  DEFENSE  PRODUCTS 

S  DYER 
85  901  AVE  53 
PO  BOX  848 
COACHELLA  CA  92236 

1  GLCC  INC 

J  RAY 

103  TRADE  ZONE  DR  STE  26C 
WEST  COLUMBIA  SC  29170 

3  PACIFIC  NORTHWEST  LAB 
M  SMITH 

G  VAN  ARSDALE 
R  SHEPPELL 
PO  BOX  999 
RICHLAND  WA  99352 


2  AMOCO  PERFORMANCE 
PRODUCTS 

M  MICHNO  JR 
J  BANISAUKAS 
4500  MCGINNIS  FERRY  RD 
ALPHARETTA  GA  30202-3944 

1  SAIC 

M  PALMER 

1410  SPRING  HILL  RD  STE  400 
MS  SH4 5 

MCLEAN  VA  22102 

1  SAIC 

G  CHRYSSOMALLIS 
3800  W  80TH  ST  STE  1090 
BLOOMINGTON  MN  55431 

1  AAI  CORPORATION 

T  G  STASTNY 
POBOX  126 

HUNT  VALLEY  MD  21030-0126 

1  APPLIED  COMPOSITES 

W  GRISCH 

333  NORTH  SIXTH  ST 
ST  CHARLES  IL  60174 

3  ALLIANT  TECHSYSTEMS  INC 
J  CONDON 

ELYNAM 

J GERHARD 

WV01  16  STATE  RT  956 

POBOX  210 

ROCKET  CENTER  WV 

26726-0210 

1  CUSTOM  ANALYTICAL 

ENG  SYS  INC 
A  ALEXANDER 
13000  TENSOR  LANE  NE 
FLINTSTONE  MD  21530 

1  OFC  DEPUTY  UNDER  SEC  DEFNS 

JAMES  THOMPSON 
1745  JEFFERSON  DAVIS  HWY 
CRYSTAL  SQ  4  STE  501 
ARLINGTON  VA  22202 
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8  ALLIANT  TECHSYSTEMS  INC 
C  CANDLAND  MN1 1  2830 
C  AAKHUS  MN11  2830 
B  SEE  MN1 1  2439 
N  VLAHAKUS  MN11  2145 
RDOHRNMN11  2830 
S  HAGLUND  MN 1 1  2439 
M  HISSONG  MN1 1  2830 
D  KAMDAR  MN1 1  2830 
600  SECOND  ST  NE 
HOPKINS  MN  55343-8367 

1  PROJECTILE  TECHNOLOGY  INC 

515  GILES  ST 

HAVRE  DE  GRACE  MD  21078 

1  LORAL  VOUGHT  SYSTEMS 

K  COOK 

1701  W  MARSHALL  DR 
GRAND  PRAIRIE  TX  7505 1 

5  AEROJET  GEN  CORP 

D  PILLASCH 
T  COULTER 
C FLYNN 
D  RUBAREZUL 
M  GREINER 

1100  WEST  HOLLYVALE  ST 
AZUSA  CA  91702-0296 

3  HEXCEL  INC 

RBOE 

F  POLICELLI 
J  POESCH 
PO  BOX  98 
MAGNA  UT  84044 

1  HERCULES  INC 

HERCULES  PLAZA 
WILMINGTON  DE  19894 

1  BRIGS  COMPANY 

J  BACKOFEN 
2668  PETERBOROUGH  ST 
HERNDON  VA  22071-2443 


1  ZERNOW  TECHNICAL  SERVICES 
L  ZERNOW 

425  W  BONITA  AVE  STE  208 
SAN  DIMAS  CA  91773 

2  OLIN  CORPORATION 
FLINCHBAUGH  DIV 
E  STEINER 

B  STEWART 
PO  BOX  127 
RED  LION  PA  17356 

1  OLIN  CORPORATION 

L  WHITMORE 
10101  NINTH  ST  NORTH 
ST  PETERSBURG  FL  33702 

1  GKN  AEROSPACE 

DOLDS 

15  STERLING  DR 
WALLINGFORD  CT  06492 

5  SIKORSKY  AIRCRAFT 

G  JACARUSO 
T  CARSTENSAN 
B  KAY 

S  GARBO  MS  S330A 
J  ADELMANN 
6900  MAIN  ST 
PO  BOX  9729 

STRATFORD  CT  06497-9729 

1  PRATT  &  WHITNEY 

C  WATSON 

400  MAIN  ST  MS  114  37 
EAST  HARTFORD  CT  06108 

1  AEROSPACE  CORP 
G  HAWKINS  M4  945 

2350  E  EL  SEGUNDO  BLVD 
EL  SEGUNDO  CA  90245 

2  CYTEC  FIBERITE 
M  LIN 

W  WEB 

1440  N  KRAEMER  BLVD 
ANAHEIM  CA  92806 
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1 

2 

2 

1 


1 

1 

1 

1 


ORGANIZATION 
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COPIES  ORGANIZATION 


HEXCEL 

TBITZER 

11711  DUBLIN  BLVD 
DUBLIN  CA  94568 

BOEING 
R  BOHLMANN 
PO  BOX  516  MC  5021322 
ST  LOUIS  MO  63166-0516 

BOEING  DFNSE  &  SPACE  GP 
W  HAMMOND  S  4X55 
J  RUSSELL  S  4X55 
PO  BOX  3707 
SEATTLE  WA  98124-2207 

BOEING  ROTORCRAFT 
P  MINGURT 
P HANDEL 
800  B  PUTNAM  BLVD 
WALLINGFORD  PA  19086 

BOEING 

DOUGLAS  PRODUCTS  DIV 
L  J  HART  SMITH 
3855  LAKEWOOD  BLVD 
D800  0019 

LONG  BEACH  CA  90846-0001 

LOCKHEED  MARTIN 
S  REEVE 
8650  COBB  DR 
D  73  62  MZ  0648 
MARIETTA  GA  30063-0648 

LOCKHEED  MARTIN 
SKUNK  WORKS 
D  FORTNEY 
1011  LOCKHEED  WAY 
PALMDALE  CA  93599-2502 

LOCKHEED  MARTIN 
R  FIELDS 
1195  IRWIN  CT 
WINTER  SPRINGS  FL  32708 

MATERIALS  SCIENCES  CORP 
B  W  ROSEN 

500  OFC  CENTER  DR  STE  250 
FT  WASHINGTON  PA  19034 


1  NORTHRUP  GRUMMAN  CORP 
ELECTRONIC  SENSORS 

&  SYSTEMS  DIV 
E  SCHOCH  MS  V  16 
1745AW  NURSERY  RD 
LINTHICUM  MD  21090 

2  NORTHROP  GRUMMAN 
ENVIRONMENTAL  PROGRAMS 
R  OSTERMAN 

A  YEN 

8900  E  WASHINGTON  BLVD 
PICO  RIVERA  CA  90660 

1  UDLP 

D  MARTIN 
PO  BOX  359 

SANTA  CLARA  CA  95052 

1  UDLP 

G  THOMAS 

PO  BOX  58123 

SANTA  CLARA  CA  95052 

2  UDLP 

R  BARRETT  MAIL  DROP  M53 
V  HORVATICH  MAIL  DROP  M53 
328  W  BROKAW  RD 
SANTA  CLARA  CA  95052-0359 

3  UDLP 

GROUND  SYSTEMS  DIVISION 
M  PEDRAZZI  MAIL  DROP  N09 
A  LEE  MAIL  DROP  Nil 
M  MACLEAN  MAIL  DROP  N06 
1205  COLEMAN  AVE 
SANTA  CLARA  CA  95052 

4  UDLP 

R  BRYNSVOLD 
PJANKEMS  170 
4800  EAST  RIVER  RD 
MINNEAPOLIS  MN  55421-1498 

1  GDLS  DIVISION 

DBARTLE 
PO  BOX  1901 
WARREN  MI  48090 
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COPIES  ORGANIZATION 


2  GDLS 

DREES 
M  PASIK 
PO  BOX  2074 
WARREN  MI  48090-2074 

1  GDLS 

MUSKEGON  OPERATIONS 
W  SOMMERS  JR 
76  GETTY  ST 
MUSKEGON  MI  49442 

1  GENERAL  DYNAMICS 
AMPHIBIOUS  SYS 
SURVIVABILITY  LEAD 
G  WALKER 

991  ANNAPOLIS  WAY 
WOODBRIDGE  VA  22191 

5  INST  FOR  ADVANCED 

TECH 
TKIEHNE 
H  FAIR 
P  SULLIVAN 
W  REINECKE 
IMCNAB 

4030  2  W  BRAKER  LN 
AUSTIN  TX  78759 

2  CIVIL  ENGR  RSCH  FOUNDATION 
PRESIDENT 

H  BERNSTEIN 
R BELLE 

1015  15THSTNWSTE600 
WASHINGTON  DC  20005 

1  ARROW  TECH  ASSO 

1233  SHELBURNE  RD  STE  D  8 
SOUTH  BURLINGTON  VT 
05403-7700 

1  R  EICHELBERGER 

CONSULTANT 
409  W  CATHERINE  ST 
BEL  AIR  MD  21014-3613 

1  UCLA  MANE  DEPT  ENGR  IV 

H  T  HAHN 

LOS  ANGELES  CA  90024-1597 


2  UNIV  OF  DAYTON 

RESEARCH  INST 
RYKIM 
AKROY 

300  COLLEGE  PARK  AVE 
DAYTON  OH  45469-0168 

1  MIT 

P  LAGACE 
77  MASS  AVE 
CAMBRIDGE  MA  01887 

1  IIT  RESEARCH  CENTER 

DROSE 
201  MILL  ST 
ROME  NY  13440-6916 

1  GA  TECH  RSCH  INST 

GA  INST  OF  TCHNLGY 
P  FRIEDERICH 
ATLANTA  GA  30392 

1  MICHIGAN  ST  UNIV 

MSM  DEPT 
RAVERILL 
3515  EB 

EAST  LANSING  MI  48824-1226 

1  UNIV  OF  KENTUCKY 

L  PENN 

763  ANDERSON  HALL 
LEXINGTON  KY  40506-0046 

1  UNIV  OF  WYOMING 

D  ADAMS 
PO  BOX  3295 
LARAMIE  WY  82071 

1  UNIV  OF  UTAH 
DEPT  OF  MECH  & 

INDUSTRIAL  ENGR 
S  SWANSON 

SALT  LAKE  CITY  UT  84 1 1 2 

2  PENN  STATE  UNIV 
RMCNITT 
CBAKIS 

212  EARTH  ENGR  SCIENCES  BLDG 
UNIVERSITY  PARK  PA  16802 
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1  PENN  STATE  UNIV 

R  S  ENGEL 

245  HAMMOND  BLDG 
UNIVERSITY  PARK  PA  16801 

1  PURDUE  UNIV 

SCHOOL  OF  AERO  &  ASTRO 
CT  SUN 

W  LAFAYETTE  IN  47907-1282 

1  STANFORD  UNIV 

DEPT  OF  AERONAUTICS 
&  AEROBALLISTICS 
S  TSAI 

DURANT  BLDG 
STANFORD  CA  94305 

1  UNIV  OF  DAYTON 

J  M  WHITNEY 
COLLEGE  PARK  AVE 
DAYTON  OH  45469-0240 

7  UNIV  OF  DELAWARE 

CTR  FOR  COMPOSITE  MTRLS 
J  GILLESPIE 
M  SANT  ARE 
G  PALMESE 
S  YARLAGADDA 
SADVANI 
D  HEIDER 
D  KUKICH 

201  SPENCER  LABORATORY 
NEWARK  DE  19716 

1  DEPT.  OF  MATERIALS 

SCIENCE  &  ENGINEERING 
UNIVERSITY  OF  ILLINOIS 
AT  URBANA  CHAMPAIGN 
I  ECONOMY 
1304  WEST  GREEN  ST 
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